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NOMENCLATURE 


A  area  (in^) 

D  forebody  diameter  (in.) 

K  specific  heat  ratio  (-) 

L  boattail  length  (in.) 

M  Mach  number  (-) 

P  pressure  ( psi ) 

r  radius  (in.) 

R  radius  of  curvature  (in.) 

Re  Reynolds  number  (-) 

5  separation  distance  measured  from  end  of  boattail  (in.) 

T  absolute  temperature  (°R) 

U*  measure  of  the  local  rate  of  acceleration  near  the  nozzle  lip 

1  [12] 

x  intrinsic  coordinate  system  (in.) 

a  angle-of-attack  (deg) 

6  boattail  angle  (deg) 

6  control  surface  deflection  angle  (deg)  or 

velocity  boundary  layer  thickness  (in.) 

<$*  displacement  thickness  of  the  boundary  layer  (in.) 

6**  momentum  thickness  of  the  boundary  layer  (in.) 

Y  specific  heat  ratio  (-) 

p  density  (lbm/in^) 

9  conical  divergence  angle  (deg)  or  momentum  thickness  of  the 
boundary  layer  (in.) 

ui  Prandtl -Meyer  angle  (deg) 


viii 


Subscripts 

b  back  or  base  conditions 

e  exit  nozzle  plane 

E  external 

F  free  jet  surface  conditions 

L  conditions  at  nozzle  lip 

m  model 

o  stagnation  state 

p  prototype 

x  flow  conditions  before  normal  shock 

y  flow  conditions  after  normal  shock 

*  throat  conditions 

«  free  stream  conditions 

*See  APPENDIX  for  nomenclature  for  a  particular  program. 
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1.  INTRODUCTION 


High  acceleration  rates  required  for  tactical  missiles  and  the  re¬ 
sulting  jet-slipstream  interaction  can  give  rise  to  undesirable  aero¬ 
dynamic  performance.  Degradation  in  performance  is  due  to  introducing 
severe  drag  penalties  through  lower  than  ambient  base  pressures  or,  as 
the  ratio  of  jet  stagnation-to-ambient  pressure  increases,  by  leading  to 
plume  induced  separation  [l,2]t,  see  Fig.  1.  Separation,  in  addition  to 
increasing  drag,  can  adversely  affect  missile  stability  and  control 
surface  effectiveness  [3], 

The  flow  component  method  of  Korst,  et  al .  [4],  has  led  to  a  basic 
understanding  of  the  problem,  including  the  establishment  of  quanti¬ 
tative  relations  accounting  for  the  influence  of  all  the  pertinent  vari¬ 
ables  and  parameters.  But  confidence  in  the  method  to  predict  actual 
separation  locations,  even  for  the  simplest  of  vehicle  geometries  [5], 
leaves  much  to  be  desired. 

For  actual  missile  or  vehicle  configurations,  especially  when  large 
angles-of-attack,  canted  fins  and  mixed  subsonic,  transonic  and  super¬ 
sonic  regions  are  present,  the  resulting  flowfields  are  very  complex, 
see  Fig.  2.  The  usual  support  of  predictive  and  corrective  prototype 
evaluations  will,  therefore,  be  strongly  dependent  on  wind  tunnel 
modeling. 

Even  then  difficult  problems  arise  in  dealing  with  the  aerodynamic 
and  thermodynamic  effects  of  propulsive  jet  plumes.  Because  it  is  in¬ 
convenient  and  often  impossible  to  carry  out  wind  tunnel  experiments 


tNumbers  in  brackets  refer  to  entries  in  the  BIBLIOGRAPHY. 


Flow  Configuration  for  Plume  Induced  Separation  from 
a  Conical  Afterbody  and  Identification  of  Geometric 
and  Operational  Parameters 


with  the  actual  hot  propellants,  one  must  revert  to  appropriate  methods 
of  plume  simulation. 

Jet  plumes  have  been  generated  in  wind  tunnels  using  a  variety  of 
methods  including  the  use  of  cold  or  heated  air  through  geometrically 
modeled  nozzles,  small  rocket  motors,  radial  injection  and  solid  bodies 
which  simulate  the  plume  shape.  Deficiencies  inherent  in  these  methods 
can  be  traced  to  failure  to  account  for  such  factors  as  specific  heat 
ratio  influence,  mass  entrainment,  viscous  effects,  wake  closure, 
temperature  and  plume  deflection.  Accounting  for  all  of  the  pertinent 
parameters  in  wind  tunnel  tests,  simultaneously,  is  impossible. 

Requirements  for  modeling  interactions  include,  in  particular: 

1.  Geometrically  congruent  inviscid  jet  contours, 

2.  Correct  pressure  rise-jet  boundary  deflection  (plume  stiffness) 
and 

3.  Mass  entrainment  characteristics  along  the  wake  boundary. 

Realizing  that  the  predominant  dynamic  interactions  depend  strongly 

on  the  specific  heat  ratio  of  the  propellant  and  to  a  lesser  degree  on 
its  temperature,  several  investigators  [6,7]  have  performed  cold  gas 
tests  using  gases  with  specific  heat  ratios  similiar  to  those  of  the 
actual  propellant  at  its  elevated  temperature.  However,  the  large 
amounts  required  and  the  relative  high  costs  of  running  with  such  gases 
makes  this  type  of  testing  often  impractical. 

A  second  approach  to  the  problem  is  to  design  test  nozzles  so  that 
the  dynamics  of  the  actual  plume  can  be  simulated  when  the  nozzles  are 
run  with  unheated  air.  Investigators  [1,8,9]  taking  this  approach  have 
used  nozzles  which  meet  the  simple  simulation  requirements  proposed  by 


5 


Goethert  and  Barnes  [10].  These  include  duplicating  the  initial  plume 
angle  as  the  jet  exits  from  the  nozzle  and  matching  the  ratio  of  jet-to- 
ambient  pressure  between  the  model  and  the  actual  (prototype)  flow. 

However,  Love  and  Lee  [11]  have  shown  that  duplicating  the  initial 
plume  angle  is  not  a  stringent  enough  requirement  to  cope  with  plume  in¬ 
duced  separation  where  the  curvature  of  the  jet  boundary  is  of  im¬ 
portance.  An  analysis  of  inviscid  jet  boundaries  near  centers  of  ex¬ 
pansion  by  Johannesen  and  Meyer  [12]  provided  a  means  for  introducing 
the  initial  plume  radius  of  curvature.  Korst  [13]  then  interpreted 
their  work  as  a  device  for  geometric  plume  modeling. 

The  assumption  of  conical  source  flow  as  an  approach  condition  to 
the  nozzle  exit,  while  greatly  simplifying  the  analysis,  becomes  in¬ 
accurate  for  nozzles  where  the  transonic  throat  region  produces  nozzle 
flowfields  which  are  significantly  different  from  those  of  ideal  conical 
source  flow.  Thus,  the  supersonic  flowfield  in  axi symmetric  nozzles  as 
affected  by  the  throat  wall  curvature  requires  careful  consideration 
[14,15,16]. 

In  addition,  replacing  a  prototype  propellant  by  an  (inert)  gas  of 
higher  specific  heat  ratio,  will  generally  lead  to  model  nozzles  with 
smaller  exit  Mach  numbers  and  exit  wall  angles.  This,  in  effect,  im¬ 
poses  a  lower  limit  on  the  range  of  geometrical  modeling  [17],  see  Fig. 
3.  The  main  problems  are  a  relatively  thick  boundary  layer  [18,19,20], 
shock  formations  [21],  and  manufacturing  nozzle  contours  with  small 
conical  divergence  angles. 

In  order  to  avoid  these  effects  and  to  extend  the  useful  range  of 
model  testing,  the  possibility  of  utilizing  model  nozzles  with  con- 


verging-djverging-converging  (CPC)  cross  sections  (including  the  con 
ventional  transonic  near  throat  geometry),  see  Fig.  4,  is  here  investi 
gated  by  analytical  and  experimental  methods. 


R/r.*8 


Figure  4a  CDC  Model  (Nozzle  A) 
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Figure  4b  CDC  Model  (Nozzle  B) 


2.  PLUME  SIMULATION  REQUIREMENTS 


Plume  simulation  is  the  attempt  to  achieve  geometrically  congruent 
jet  contours  between  the  model  and  prototype  flows  with  the  most  perti¬ 
nent  dynamic  characteri sties  satisfied. 

2.1  Geometric  Simulation 

Previous  techniques  have  been  inadequate  in  that  they  only  du¬ 
plicate  the  initial  plume  angle.  The  current  method  matches,  in  ad¬ 
dition,  the  initial  radius  of  curvature  needed  to  adequately  define  the 
plume  shape  bordering  the  near  wake.  With  reference  to  Fig.  5,  this  can 
be  expressed  as 


(1.2) 


where  m  denotes  the  model  flow  and  p  the  prototype  flow.  It  is  assumed 
that  the  following  parameters  are  given:  Ypp,  MLp,  9Lp  and  Mpp. 

The  geometric  simulation  parameters,  0p  and  Rp,  are  related  to 
flow  conditions  at  the  nozzle  lip  by  a  second-order  series  solution  of 
the  axisymmetric  compressible  potential  flow  equations  near  centers  of 
expansion  [12].  As  a  result  of  such  an  analysis,  a  model  nozzle 
configuration  can  be  determined  (for  either  ideal  or  non-ideal  conical 
source  flow)  [17].  It  is  important  to  note  that  geometrical  plume 
modeling  accomplished  by  this  analysis,  while  establishing  proper  values 
for  initial  expansion  angle  and  plume  radius  of  curvature,  can  be  inter¬ 
preted  to  be  a  second  order  initial  condition  for  a  more  accurate  deter¬ 
mination  of  plume  contours  by  the  Method  of  Characteristics  (MOC).  In 


addition,  it  is  significant  that  satisfying  the  geometrical  congruence 
of  the  plume  contour  still  leaves  one  free  parameter  available  for 
achieving  dynamic  simulation. 

2.2  Dynamic  Simulation 

Ideally  the  dynamic  specifying  condition  should  account  for  both 
inviscid  and  viscid  aspects  of  the  plume.  This  apparently  requires  two 
closure  conditions  to  be  sati sfied--the  recompression  ratio  at  the  end 
of  the  wake  and  conservation  of  mass  within  the  wake.  Due  to  geo¬ 
metrical  requirements,  only  one  parameter  is,  however,  available. 

For  studying  plume  induced  separation,  it  is  recommended  [22]  that 
plume  surface  pliability  be  modeled  to  satisfy  both  the  inviscid  and 
viscid  interaction  mechanisms  of  plume-slipstream  confluence  at  the  end 
of  the  wake.  The  selection  of  inviscid  pressure  rise-deflection  re¬ 
lations  along  the  plume  boundary  was  chosen  because  the  jet  plume  es¬ 
sentially  contributes  to  the  wake  closure  through  its  inviscid  jet 
boundary  geometries  before  and  after  impingement.  The  viscid  aspects  of 
jet  mixing,  attributed  predominately  to  the  separated  slipstream,  are 
diminished  by  the  relaminarization  of  the  nozzle  boundary  layer  during 
its  rapid  expansion  after  the  nozzle  exit. 

For  relatively  small  angular  deflections  (beyond  the  linearization 
limits),  the  requirement  for  plume  pliability  (jet  surface  Mach  number) 
is  found  from  Weak  Shock  Modeling  (Program  WSHOCK)*. 


tProgram  listings  are  found  in  the  APPENDIX. 
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2.3  Calculation  Procedure 

One  intends  to  model  a  prototype  plume  that  results  from  y 
MLp,9Lp  and  MFp*  yie1din9  an  initial  slope  9pp  and  an  initial  radius  of 
curvature  Rpp.  For  a  given  model  gas  vm»  using  either  Weak  or  Strong 
Shock  Modeling  will  determine  Mpm. 

For  a  selected  trial  value  of  M|_m,  the  Prandtl -Meyer  relation 
for  9^  and  the  analysis  of  Johanessen  and  Meyer  [12]  will  yield  the 
radius  of  curvature  ^Fm*  If  a  solution  exists,  it  will  be  found  by 
satisfying  Rpm  =  RFp. 

Nozzle  geometries,  nozzle  Mach  numbers,  and  jet-to-ambient  pressure 
ratios  can  differ  widely  from  model  and  prototype.  These  parametric 
variations  do  not  seem  to  impair  the  accuracy  of  matching  the  plume 
surface  geometry. 

Due  to  the  complexity  of  the  iterative  calculation  procedures,  one 
must  revert  to  computational  solutions. 


3.  FLOW  MODEL 


Although  geometric  simulation  requirements  involve  only  inviscid 
flow  mechanisms,  the  inherent  dependence  of  the  problem  on  viscous  mass 
interchange  requires  viscid-inviscid  analysis  for  adequate  closure. 

3.1  Nozzle  Flow 

3.1.1  Inviscid  Flowfield 

The  inviscid  analysis  of  the  flow  in  the  nozzle  consists  of 
a  series  solution  of  the  transonic  flow  in  the  vicinity  of  the  throat 
and  an  inviscid  MOC  solution  for  the  supersonic  flow  downstream  of  the 
throat  region. 

The  disturbing  effects  of  throat  wall  curvature  on  the  transonic 
flow  region  require  an  accurate  initial  condition  to  properly  calculate 
the  supersonic  flowfield.  Transonic  flow  analyses  by  Sauer  [15], 
Kliegel  and  Levine  [16],  or  Dutton  and  Addy  [23]  can  be  used  to  provide 
the  noncharacteristic  initial  condition  for  the  MOC  nozzle  solution. 

The  MOC  (Program  CDCSECT)  is  then  used  to  calculate  the  supersonic 
flowfield.  Nozzle  wall  geometry  and  model  propellant  properties  serve 
as  input.  Characteristics  are  calculated  throughout  the  nozzle's  region 
of  influence  on  the  developing  plume. 

The  exit  conditions  of  the  inviscid  nozzle  flow  serve  as  input  to 
the  MOC  plume  solution,  see  Section  3.2.1,  and  are  also  used  in  deter¬ 
mining  the  second-order  modeling  scheme  for  non-ideal  conical  source 


3.1.2  Viscous  Boundary  Layer  Effects 

The  viscid  analysis  of  the  nozzle  flow  is  an  extension  of 
established  methods  for  compressible  laminar-transitional-turbulent 
boundary  layer  development  [19,24]  (Program  GRTRAN1). 

Beginning  with  the  Gruschwitz  [25]  method  for  a  laminar  com¬ 
pressible  boundary  layer,  the  transition  to  turbulent  flow  is  based  on  a 
criterion  related  to  the  local  Reynolds  number  based  on  an  equivalent 
displacement  thickness  as  a  function  of  the  local  second-shape  factor 
[19].  The  turbulent  portion  of  the  boundary  layer  was  then  calculated 
using  a  scheme  based  on  the  method  of  Culick  and  Hill  [24].  No  dis¬ 
tinction  has  here  been  made  between  the  point  of  instability  and  the 
actual  onset  of  turbulence. 

Since  it  was  found  that  the  nozzle  boundary  layer  was  predominately 
turbulent  (transition  occurred  before  the  throat),  comparison  with  the 
fully  turbulent  results  based  on  the  Bartz  method  [18]  (Program  BARTZ) 
was  attractive.  Close  agreement  between  the  two  can  be  observed,  see 
Fig.  6. 

Both  methods  utilize  the  approximate  integral  momentum  and  energy 
boundary  layer  equations  which  are  solved  for  local  velocity,  displace¬ 
ment,  momentum  and  energy  thicknesses.  As  part  of  the  output,  the  BARTZ 
program  calculates  a  Reynolds  number  (based  on  the  velocity  boundary 
layer  thickness)  that  is  used  in  the  calculation  of  the  free  shear  layer 
development,  see  Section  3.2.2,  and  in  the  calculation  of  plume  impinge¬ 
ment,  see  Section  5.3. 

The  boundary  layer  displacement  thicknesses  were  small  enough  so 
that  an  iterative  procedure  for  obtaining  improved  inviscid  nozzle  flow 
field  solutions  was  not  deemed  necessary. 
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3.2  Plume  Development 

3.2.1  Inviscid  Plume  Boundary 

The  MOC  (Program  JETC)  is  used  to  calculate  the  expanded  in- 
viscid  plume  boundary.  The  final  characteristic  line  calculated  by  the 
MOC  nozzle  solution,  see  Section  3.1.1,  is  used  as  a  starting  condition 
for  the  expansion  centered  at  the  nozzle  exit.  Jet  surface  Mach  number 
is  specified  to  account  for  plume  pliability,  see  Section  2.2.  Circular 
arcs  are  then  fitted  between  the  characteristic  boundary  points  to  ac¬ 
count  for  angular  variation  in  the  plume  surface. 

A  second-order  modeling  scheme  for  non-ideal  conical  source  flow 
(Program  NOCON)  is  used  as  a  comparison,  see  Fig.  7.  Flow  parameters  at 
the  nozzle  exit  (ML,  eL,  U^)  calculated  by  the  MOC  nozzle  solution, 
serve  as  input. 

3.2.2  Viscous  Mass  Entrainment 

While  it  can  be  expected  in  cases  involving  large  regions 
of  plume  induced  separation  that  the  mass  and  energy  balance  in  the  wake 
is  predominantly  controlled  by  the  viscous  mechanisms  of  the  slipstream 
(external  flow),  it  is  nevertheless  of  interest  to  investigate  the 
viscous  mechanisms  along  the  jet  boundary.  In  absence  of  a  confluent 
slipstream,  plume  impingement  of  solid  walls  will  indeed  be  entirely 
controlled  by  the  vi scid-invi scid  interactions  of  the  plume  alone 
[26],  Mass  bleed  may  be  necessary  for  the  proper  modeling  of  such  im¬ 
pingement  flows  [27,28]. 


A  computer  program  (Program  FREES)  is  available  for  determining  the 
development  of  the  free  shear  layer  in  the  form  of  a  two-layer  com¬ 
pressible  laminar-transitional-turbulent  model  based  on  local  similarity 
concepts  [29],  see  Fig.  8.  While  it  is  not  possible  to  make  accurate 
measurements  from  the  Schlieren  photograghs.  Figs.  9,  there  is  general 
agreement  between  theory  and  experiments. 
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Figure  8  Free  Shear  Layer  Development 
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4.  EXPERIMENTAL  INVESTIGATION 


4.1  Objectives 

The  objective  of  this  thesis  is  to  determine  the  feasibility  of  de¬ 
signing  and  operating  a  CDC  nozzle  to  extend  the  range  of  plume  simu¬ 
lation.  To  do  so,  the  following  were  designated  as  specific  objectives: 

1.  Establish  limits  for  conical -to-conical  nozzle  modeling, 

2.  Design  a  nozzle  based  on  theoretical  analysis  and  explore  the 
extension  of  these  limits  by  use  of  CDC  nozzles, 

3.  Experimentally  verify  the  theoretical  calculations  including 
starting  characteristics  and  flowfield  peculiarities  (coa¬ 
lescence  of  disturbance  lines,  etc.),  and 

4.  Determine  the  potential  and  limitations  of  using  such  nozzles 
as  modeling  tools. 

4.2  Apparatus  and  Instrumentation 

Two  specific  nozzle  configurations  were  selected  from  a  variety  of 
candidate  configurations.  The  reason  for  the  selection  was  that  the  two 
produced  the  intended  plume  shape  and  had  an  almost  constant  Mach  number 
wall  distribution  after  the  conical  supersonic  expansion  section.  This 
allowed  changing  the  exit  lip  angle,  9^,  without  a  change  in  exit  lip 
Mach  number.  The  CDO  configuration  was  derived  from  the  idea  of 
replacing  the  conventional  conical  supersonic  expansion  section  and  free 
jet  boundary  by  an  appropriate  combination  of  conical  and  circular  arc 
wall  sections. 

The  two  nozzles,  see  Fig.  4,  were  fabricated  to  study  the  effect  of 
throat  wall  curvature.  This  effect  is  evident  in  Fig.  10.  As  seen  in 


Figs.  11  and  12,  application  of  the  Kliegel  and  Levine  [16]  throat  so¬ 
lution  successfully  predicts  the  transonic  flowfield. 

Since  the  analytical  design  process  was  done  non-dimensional ly,  the 
size  of  the  throat  determined  all  other  dimensions.  An  analysis  of  the 
compressor/storage  system  led  to  a  throat  radius  of  0.375  inches. 

The  wall  pressure  data  and  Schlieren  photographs  were  taken  with 
the  nozzles  exhausting  into  the  ambient.  Pressure  measurements  were 
taken  with  a  conventional  transducer-scani val ve-ampl i fier  setup.  An 
Apple  Il-e  computer  was  employed  as  a  data  reduction/recording  device. 

For  the  impingement  and  stagnation  pressure  data,  an  enclosure  was 
fabricated.  The  enclosure  served  several  uses,  it  positioned  the  ex¬ 
haust  tube  that  the  impingement  data  was  taken  from,  it  secured  the 
pitot  probe,  it  allowed  evacuation  of  the  test  section  (lowering  the 
back  pressure)  and  lastly  it  significantly  decreased  the  noise  level 
during  tests.  If  mass  bleed  was  necessary,  a  metering  valve  could 
easily  be  connected  to  the  enclosure. 

Pressure  taps  in  the  exhaust  tube  provided  a  means  for  collecting 
impingement  data.  The  inner  radius  of  the  tube  was  1.4  times  larger 
than  the  nozzle  exit  radius.  This  allowed  plume  expansion  with  the 
assurance  of  impingement. 

All  testing  was  done  at  the  Gas  Dynamics  Laboratory,  Department  of 
Mechanical  and  Industrial  Engineering,  University  of  Illinois  at  llrbana- 
Champaign. 


?>v>v 


M.O.C.  (Program 
CDCSECT) 
Experimental  Data 



Distribution 


26 


J 


5.  EXPERIMENTAL  RESULTS 
5.1  Nozzle  Wall  Pressure  Distribution 


Nozzle  wall  pressure  data  taken  from  the  models  are  shown  in 


Figs.  11  and  12  compared  to  the  MOC  (Program  CDCSECT)  and  one¬ 


dimensional  solutions.  Note  the  constancy  of  the  Mach  number  after  the 
conical  divergence  section. 

5.2  Plume  Development 

5.2.1  Plume  Shape  Comparison 

The  Schlieren  photograph.  Fig.  9c,  was  used  to  compare  the 
inviscid  solutions  of  the  MOC  (Program  JETC)  and  of  the  second-order, 
non-ideal  conical  source  flow  approximation  (Program  NOCON).  The  com¬ 
parison  is  shown  in  Fig.  7. 

In  Fig.  9,  note  the  disturbance  lines  emerging  from  inside  the 
nozzle.  These  were  predicted  and  can  be  shown  graphically  on  a  plotting 
version  of  the  MOC  nozzle  solution  (Program  CDCSECT).  The  disturbance 
lines  are  focused  by  the  concave  wall  curvature  of  the  nozzle  and 
eventually  coalesce  into  oblique  shock  waves.  Neither  program,  JETC  or 
NOCON,  takes  these  disturbance  lines  into  account  in  calculating  plume 
shapes.  Even  so,  agreement  between  calculated  and  observed  plume 
boundaries  is  satisfactory,  see  Fig.  7. 

Even  though  the  coalescence  of  compressive  disturbance  lines  cause 
discontinuities  in  plume  slope,  this  is  detrimental  only  if  the  coa¬ 
lescence  occurs  upstream  of  the  confluence  with  the  external  slipstream 
or  if  it  interferes  with  the  plume  boundary  in  the  wake  closure  region. 
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5.2.2  Mach  Number  Survey 

The  survey  was  done  by  measuring  the  Rayleigh-pitot  stag¬ 
nation  pressure  using  a  0.050  in.  diameter  pitot  probe  and  using  normal 
shock  relations  to  determine  the  Mach  number  Mx.  Figure  13  shows  the 
comparison  of  experimental  data  with  the  inviscid  solutions  (Programs 
COCSECT  and  SECTP). 

Non-parallel  outflow,  shock  wave  interference,  and  mechanical  in¬ 
accuracy  are  some  of  the  reasons  why  the  theory  does  not  better  predict 
the  Mach  number  distribution. 
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5.3  Jet  Impingement  on  Solid  Cylindrical  Surface 

The  impingement  study  was  done  to  further  characterize  the  modeling 
capabilities  (or  limitations)  of  the  CDC  configurations.  The  theo¬ 
retical  calculations  (Program  ASUEXP)  have  been  previously  determined  to 
be  reliable  [30]  for  predicting  viscid  and  inviscid  impingement  aspects 
of  conventional  converging-diverging  nozzles  discharging  into  sudden  en¬ 
largements  in  cross  section. 

In  the  present  study,  predicted  and  measured  pressure  distri¬ 
butions,  see  Fig.  14,  were  analyzed.  While  in  agreement  for  base 
pressure,  they  exhibit  distinctive  differences  near  and  after  plume  im¬ 
pingement.  This  may  be  attributable  to  the  proximity  of  the  coalescent 
compression  waves,  which  are  presently  unaccounted  for  by  the  inviscid 
jet  plume  boundary  calculations.  The  calculated  length  of  the  viscid 
interaction  region,  (Xinviscid"xviscid^re»  can  be  considered  as  a 
measure  for  the  pressure  feedback  domain  through  the  impingement  wall 
boundary  layer.  This  serves  as  an  illustration  for  the  limitations  in¬ 
herent  in  the  use  of  CDC  nozzles  alluded  to  in  Section  5.2.1. 
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6.  CONCLUSIONS 

1.  Success  of  the  general  modeling  scheme  has  been  established  in  as 
much  as  the  method  allows  simulation  of  aerodynamic  plume  inter¬ 
action  effects  of  great  detail  even  under  complex  operational  con¬ 
ditions  (large  angles-of-attack,  boattail,  etc.),  see  Fig.  15. 

2.  The  CDC  nozzle  configurations  do  extend  the  range  of  geometrical 
modeling,  see  Fig.  16,  allowing  the  modeling  of  prototype  nozzles 
with  lower  exit  lip  Mach  numbers  using  model  propellants  with  higher 
specific  heat  ratios. 

3.  The  wall  curvature  of  the  end  section  of  the  CDC  nozzle  produces 
convergence  of  compressive  lines  which  leads  to  plume  boundary  de¬ 
flections  downstream  of  the  nozzle,  this  limits  the  modeling  range 
unless  large  radii  of  curvature,  R/r*>8,  are  used  or  these  effects 
are  downstream  or  small  in  the  region  of  interest. 

4.  Agreement  between  predicted  and  experimentally  determined  flow 
fields  (inviscid  and  viscid),  within  the  limits  stated  above,  has 
been  observed. 

5.  While  plume-slipstream  interactions  have  not  been  carried  out  in 
this  investigation,  it  can  be  expected  that  the  successful  extension 
of  the  modeling  methodology  can  be  established  by  experimentation  in 
appropriate  wind  tunnel  facilities  [31], 
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Figure  15  Oil  Flow  and  Schlieren  Photographs  at  ot  =  -10° 
and  PQ  =  1.00  MPa  for  the  Freon  (Prototype)  and 

Air  (Model)  Nozzles 

(a)  Freon  M.  =  2.6 

(b)  Air  M,  =  1.41 
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APPENDIX 


This  APPENDIX  is  a  compilation  of  all  of  the  computer  programs  that 
were  used  in  this  thesis.  Almost  all  of  the  theoretical  work  done  for 
the  thesis  has  been  incorporated  into  these  programs  which  were  de¬ 
veloped  in  the  Department  of  Mechanical  and  Industrial  Engineering,  Uni¬ 
versity  of  Illinois  at  Urbana-Champaign,  Urbana,  Illinois.  Because  of 
the  number  and  complexity  of  the  programs,  each  is  listed  with  a  brief 
discussion  of  what  the  program  does,  the  input  and  output  and  pertinent 
nomenclature. 

The  listing  includes  the  following  programs: 
program  page 


CDCSECT 

34 

JETC 

46 

NOCON 

53 

WSHOCK 

56 

SECTP 

58 

PSIM 

67 

GRTRAN1 

72 

BARTZ 

80 

FREES 

86 

ASUEXP 

91 

All  programs,  unless  otherwise  stated,  are  written  in  BASIC  and 
were  run  on  the  CY8ER  175.  As  an  aid  in  understanding  what  each  program 
does  and  its  relationship  to  the  thesis  and  other  programs,  a  flowchart 
is  included  on  the  following  page. 


^BOUNDARY  LAYER 


■FREE  SHEAR 


-IMPINGEMENT-*) 


LAYER  DEVELOPMENT 


Program  CDCSECT 


Program  CDCSECT  calculates  the  inviscid  supersonic  nozzle  flowfield 
by  the  MOC  after  the  transonic  flow  region  has  been  calculated  by  a 
method  based  on  the  analysis  of  Kliegel  and  Levine  [16],  Individual 
characteristics  are  printed  and  the  last  characteristic  line  is  stored 
in  file  "SONIC"  to  be  used  as  input  for  Programs  JETC  and  SECTP.  In 
addition,  a  routine  has  been  added  to  allow  the  calculation  of  flow 
properties  at  a  given  cross  section. 


input 

quantity 

symbol 

explanation 

y 

GO 

specific  heat  ratio 

re/r* 

00 

nozzle  exit  radius* 

R*/  r* 

CO 

throat  radius  of  curvature, 

(see  Fig.  4) 

9 

TO 

conical  wall  divergence  angle  (deg) 

k 

number  of  points  in  the  MOC  net 

X4 

x  location  of  the  end  of  the  conical 

divergence  section,  measured  from 

the  throat 

R/r* 

UO 

radius  of  curvature  for  the  circular 

arc 

LO 

x  location,  from  the  nozzle  exit, 

for  flow  profile  determination 

tAl 1  nozzle  dimensions  are  non-dimensional i zed  by  the  throat  radius. 


70  FILE  Ml  ="S0NIC“ 

(30  Dirt  A<(3),B4(3),C*(3> 

90  BI="Y£S" 

100  GOTO  5140 

110  IF  X 3 > E  AND  P<0  THEN  2130 
I  20  Y3=X2 

130  Y4=X2  +  ( UO* ( SQR ( 1  -  <  ( Z2-Y3 )./U0)  2 ) -C05<  TO  >  >  +  23-R2 > /T? 

140  IF  ABS(  Y4-Y3K0.0001  THEN  170 
150  Y3=Y4 
I  60  GOTO  130 
170  X  3  =  Y  4 

180  R3=R2+ (X3-X2) *T9 
190  T3=FNZ ( ( 22-X3 ) /UO ) 

200  GOTO  2160 

210  DIM  X(9,1l >,R< 9,11 ),H< 9,11 ),T<?,11 ),DC3,14),E (3,14), F'3,1 4) 
720  DIM  G(3,35),K(7>,H<7),I(7>,J(7>,V(4,20) 

; ’ 30  DIM  P ( 1 , 1 1 , 0 f  1 , 1 1 ) , B ( * , !  *  : , C  !  1  . 1  1 ,  A ( 1  5 , !  1 ) 

‘ 40  2EF  FWA(X)  =  ATN(SQR(1-X*X)/X) 

230  DEF  FNA(X)=ATN(1/SQR(X*X-1 ) ) 

760  DEF  FNS(X>=SQR(2>*X/SQR(GO*1-'GO-1  >  «=X ' 

770  DEF  FNQ(A)*1/<M*TAN(A) ) 

280  DEF  FNF ( A) =SIN( T)*SIN( A),  SIN< T+A ) 

290  DEF  FNG(A)=SIN(T)tSIN(A)/SIN(T-A) 

300  DEF  FNZ ( X ) =ATN(X/SQR( 1  - X * X )  > 

310  PRINT  "CONICAL  NOZZLE  SECTION  STARTS  AT  X="X7;"R="R7 
320  PRINT  "FOR  CDC  NOZZLE  INPUT  X4>X7  0="; 

330  INPUT  X4 

340  IF  X 4  =  0  THEN  420 

350  Z3=R7+(X4-X7)*TAN(T0) 

360  PRINT  "ENTER  RADIUS  OF  CUR V., CDC  SECTION-"; 

170  INPUT  UO 

380  PRINT  "CURVED  WALL  STARTS  AT  X*“X4"A'  PRC IU$~"73"R 0  !US  CF 
;o /}  22- /J t-MO » :  I!J ■  T0  > 

■100  PRINT  "MAX.  RADIUS  OF  NOZZLE8 "Z3+U0* ( 1 -COS< TO) ) "AT  X-"Z2 
410  GOTO  6060 
420  I=S=P=U=1 

430  PRINT  "NOZZLE  EXIT  AT  X="E 
440  Q=N=Z0=0 
450  Z-M 

460  PRINT  "INITIAL  POINTS  ALONG  NON -CHARACTERISTIC  CURVE" 

470  PRINT  "1=1" 

480  FOR  J=1  TO  K 

490  PRINT  J,X(I,J),R(I,J),M(I,J),T(I,,J) 

300  NEXT  J 


era 


510  PRINT 
520  PRINT 

530  PRINT  “FLOUFIELD  COMPUTED,  ALL  POINTS  PRINTED" 

540  PRINT 

550  GOTO  740 

540  FOR  J=1  TO  K-S 

570  X I =  X ( I , J ) 

580  X 2  =  X ( I , J+ 1  ) 

590  R 1 =  R ( I , J ) 

600  R2=R( I, J* 1 ) 

610  Ml =M( I , J) 

620  M2=M( I, J+1 > 

630  T*  =T( I, J) 

640  T2=T(I, J*  ■ 

6Z0  jOSUB  930 
iiO  T!I+1,J>;T3 
670  ?  =p 7 

.430  M'I+1, J>=H3 
690  T(I+1,J>=T3 

695  IF  X 3 : L 0  AND  XI  LO  THEN  597 
•.'j 96  GOTO  705 
697  GOSUB  4475 

705  IF  X 3  L 0  AND  X 2  L 0  THEN  707 
704  GOTO  720 
70 7  GOSUB  6730 
720  NEXT  J 
’30  RETURN 
7 40  GOSUB  540 
750  IF  [=1  THEN  810 
740  IF  Z>0  THEN  810 
770  IF  P<0  THEN  6190 
780  GOSUB  1350 
790  GOSUB  2210 
800  GOSUB  1350 
810  GOSUB  2030 
1340  1  =  1  +  1 
1370  S=S+t 
880  U-=U+1 
890  Z=-Z 

900  IF  I <3  THEN  920 
910  GOSUB  5910 
920  GOTO  740 
930  M9=M1 

940  A1=FNA(FNS(M1 ) > 

950  A2=FNA<PFNS (H2> ) 

940  A4=A 1 
970  A 5  =  A 2 
930  T4=T1 
990  T5-T2 
1000  M4-M1 
1010  H5=N2 
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1020  R4=R1 
1030  R5=R2 

1040  X3*<R2-R1*yi*TAN(T4-A4>-V2*TAN<T'i+A5>' /( fAN<Tl-A4'  TAN(T5+A5) > 
1050  R3*R1 +(X3-X1)*TAN(T4-A4) 

1060  T=T4 
1070  A=A4 
1080  H*H4 
1040  Q1 =FNQ(A) 

1100  G 1 =FNG ( A ) 

1110  T=T5 
1120  A=A5 
1110  N=N5 
1140  Q2-:FNQ(  A ) 

1150  IF  R2*0  F HEX  2180 
1160  F 2  =  FN F (A) 

1170  M3=(<  P-T2)+Gt*(R3-R1  VR4+F2*' R3*R2)/R5*Q1  *H1  *02+«2)/(Q1  *-02) 
1180  IF  M3> 1  THEN  1200 
1190  H3=M1 

1200  T7=Q2*(M3-«2)-F2*1R3-R2)/R5 
1210  T3=T2+T7 

1220  IF  ABS(N?-N3)\0.000t  THEN  1340 

1230  H9=M3 

1240  T4=m+T3>/2 

1250  T5=(T2+T3)/2 

1260  N4»(H1+R3>/2 

1  270  N5  =  < H2+H3)  '2 

1280  A3«FNA(FNS(H3>> 

1290  A4=(A1+A3)/2 
1300  A 5= ( A2+A3 ) /2 
1310  R5=(R3+R2)/2 
1320  9 'R3+R1 1/2 
’320  GOTO  1040 
13 jo  F'ETiJRN 
.'35 0  y2-'< (1,1 ) 

1  360  R 2 = R ( 1,1) 

1370  M2=M(I,1> 

1380  T2«T(If1) 

1390  A2  =  FNA ( FNS  <  H2 ) ) 

1400  N9=«2 
1410  A5=A2 
1420  T5=T2 
1430  H5~N2 
1440  R5-:R2 
1450  T?=TAN<T5»A5) 

1460  IF  X2>X7  THEN  2120 
! 470  IF  X 3 > X 7  THEN  2120 
1480  IF  *M"YESM  THEN  1600 
I  490  X0=< <X2*T9-R2)+T?)/(T?+T9-C0> 

1500  X?=(R2*R2+X2*X2*T9*T9-2*R2*X2*T9-1 )/(T9+T9-C0> 

1510  X3=+X0+SQR(X0*X0-X9) 

1  520  IF  X 3 > X 7  THEN  2120 


1530  R3=SQR(1+C0*X3*X3) 

1540  T3=C0*X3/SQR( 1 +CO*X3*X3) 

1550  T3=ATN(T3) 

1560  Q2=FNG(A5) 

1570  F2=FNF(A5) 

1580  H3=M2+(T3-T2+F2*(R3-R2)/R5)/Q2 
1590  IF  BIO" YES"  THEN  1700 
1600  X 0  =  C 0  + 1  -R2+T9*X2 
1410  X9=T9*X0/M+T?*T9) 

1-420  X3=X9-SQR ( X9*X9+(CO*CO-XO  *:XO)/( T9+T9+1 ) ) 

1430  IF  X3/X7  THEN  2120 

1440  R3=1+CO*(1-SQR(1-X3tX3/(CO+CO)>) 

1650  T3*ATN( 1/SQR(C0*C0/(X3*X3) -1 )) 

1  660  PRINT  X3,R'3,T3*57.2958 
1670  Q2=FNQ(A5) 

1680  F2=FNF< A5) 

1  690  M3=H2+  <  T3-T2+F2*  <  R3-R2 ) /R5 )/Q2 

1700  IF  ABSIM9-H3X0.0001  THEN  1780 

1710  N9=N3 

1720  T5=(T2+T3)/2 

1730  H5*(N2+N3)/2 

1740  A3=FNA<FNS(f13)  > 

1750  A5=(A2+A3)/2 
1760  R 5= <  R 3  «-Fi2 )  / 2 
1770  GOTO  1450 
I  780  X 8 = X 3 
1790  R8=R3 
1800  H8=N3 
1810  T8=T3 

1  820  IF  X3H0  ANP  Y2  10  THEN  1322 

1821  GOTO  1840 

1822  GOSUB  6730 
1840  RETURN 

1850  FOR  -'C  t  STEP  ■'-P 

I860  V f  I  + 1  r .! ■*■  1  ■  =  *•  ,  J> 

1  3^0  R(  1  +  1 ,  J  +  1  )  =  R(IM  ,J) 

1880  H(I+t,J*1>«H<I+1,J> 

1890  T( 1  +  1 ,  J+1 1  * T ( I ♦ 1 , J ) 

1900  NEXT  J 
1910  J=1 

1920  X(I*\1)*X8 

1930  R( 1  +  1 ,1 ) =R8 

1940  H(I+1,!)»H8 

1950  T( 1*1 , 1 )=T8 

I960  IF  X(  1  +  1 ,1  XX7  THEN  1990 

1970  IF  P<0  THEN  1990 

1980  GOSUB  2570 

1990  K=K+2 

2000  IF  P<0  THEN  2020 
2010  IF  X( 1+1 , 1 )>E  THEN  2740 
2020  RETURN 


2030 

2040 

2050 

2060 

2070 

2080 

7100 

2110 

2120 

7130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

3270 

73SO 

2390 

3400 

2410 

.3420 

'7430 

2440 

2450 

2460 

2470 

2471 

?47'j 

2490 

2500 

2510 

2520 

2530 


IF  1=2  AND  P<0  THEN  3020 

PRINT  "KOU  »U="U 

FOR  J=1  TO  K-S 

IF  J<Q  THEN  2100 

IF  J-Q  AND  P  .0  THEN  3040 

PRINT  J,X(I*t,J),R(I*'1»J>»N(I  +  1»J)»T'I  +  1lJ)  *-57.2958 

NEXT  J 

RETURN 

IF  X 3 > X 4  AND  X 4 '  0  THEN  *10 

X3=  (R2-F;7  +  X7*TAN(  TO  >  -  X2  *T?)/(ThN(  TO)  -  T9) 

R3=R7+(X3-X7)*TAN(T0) 

T3  =  T0 

IF  B$="YES"  THEN  1670 
GOTO  1560 

H3*(Q2*H2/2+T1+Q1*H1  +G1  ■*  'R.3-R1  -■  ,'R4 ) / <32/2+01  1 
T7=Q2*(H3-N2)/2 
GOTO  1210 
XT  =X< I , K-S  +  1  ) 

R1 =R( I , K-S+1 ) 
ill  =ft(  I , K - S+ 1 ) 

T1 =T( I , K-S+1 ) 

A 1 =FN A ( FNS ( Ml ) ) 

N9=1 .01 

A4=A1 

T4=TI/2 

M4=M1 

R4=R1/2 

T9=TAN(A4-T4) 

X  3  =  X 1 +R1/T9 
G1=FNG<A4) 

Q1 =FNQ (A4 ) 

N3=N1+(T1-G1 )/3t 
PRINT  "h3,AXIS-“fl3 
IF  ABS<H3-M9) <0.0001  THEN  2430 
H9  =  M3 

N4  =  (M1 +  M  3 ) / 2 
A3  =  CNA 1 TN3 ( M3 )  • 

44= ''  A1  +A3)  .'2 

goto  ::r  ' 

XC 1+1 , K-S+1 )*X3 

R<  I  +  T ,  K-S  +  1  >  --0 

H(I  +  1 f  K-S  +  1 »=H3 

T(I+1 , K-S+1 )=0 

IF  X 3 :  L 0  AND  XI <L0  THEN  2472 

GOTO  2490 


RETURN 

FOR  J=l  TO  K-S 
X  < 1 , J)=X(3,J) 
R<1 ,J)=R(3,J) 
N(  1 ,  J)=(1(3f  J) 
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T(  1 , J)=T(3, J) 

NEXT  J 
GOTO  2020 

D1 =FNA(FNS(M( 1-1 , 1 ) ) ) 

D2=FNA(FNS(H( 1,1))) 

D3=-C0S(D1 ) +C0S(D2 )*H( I , 1 )/H( 1-1 f 1 ) 

D4=SQR ( £X< 1-1 ,1 ) -X  < 1 , 1 ) >*2+<R(I-1 ,1 >-R(I,l ))‘2)/R(T-1 , ! ) 

U1=D3*N(I-1,1)/B4 

PRINT  "U1*="U1 

D5=D3/D4+SIN«!  TO  >  * < 5IN( 02) )  2 
D  6 = D  5 / (2  +  (CCS ' 52 ' 1 ' 2> 

PRINT  "ACCELERATION  TEF:f1="D5 

PRINT  "AT  WALL  POINT  WHERE  R/R*»"R( 1-1,1) "AND  X/R* ="X( I- 1 , 1 ) 

PRINT 

PRINT 

RETURN 

PRINT  "EXIT  CONDITIONS  FOR  HM"FNS(H8> 

PRINT  "AT  X  =  "X8; "R="R8,* "UHERE  THETA*  "TB+'i?.  2958"  DEG" 

B6=X8 

B7=R8 

BS=M8 

B  9  =  T  8 

PRINT 

GOTO  3000 

01*0 

GOSUB  2950 

UO«(E-X(I-f , 1 ) )/(X(l+1 , ! )-X(I-1 ,  m 
PRINT  00 

HI  =X  ( 1-1  r  1  )  +  ( X ( 1 , 1 ) -X ’  I  -1 , 1 > !  HJO 
H2=R( 1-1 , 1 )+(R( I, 1 )-R( 1-1 , 1 ) ) *U0 
H3  =  N(  1-1 , 1  J  +  ( N( 1, 1 ) - M ( 1  —  1 , 1  ■  )*U0 
E4 =T  ( I  - 1 , 1 )  +  ( T  ( 1 , 1 )  -T  ( I  - 1 , 1  •  '  < 910 
X  ( 1 ,  1 >  =H1 
?.<  I,'  )  =H2 
7 :  lr  ‘ 

H(I,1 )=H3 
GOSUB  1 350 

IF  ABS-U0-U1 XO.OOOOI  THEN  2700 

R(I  +  1 ,1  )  =  R8 

x •: i  + 1  ,i  )=X9 

H(I  +  1 ,1 )=H8 

T(  1  +  1 f  1  >*T8 

U1  =U0 

GOTO  2750 

GOTO  2700 

V1=X(I,1) 

V2=R  <  1 , 1 ) 

V3=T (1,1) 

V4 =H ( 1 , 1 ) 

A6=V1 

A7=V2 


45 


2983  A8=V3 

2984  A9=V4 
2990  RETURN 
3000  P=-1 
3010  GOTO  2020 
3020  Q=Q+ 1 
3030  GOTO  2040 
3040  N=N  +  1 

3050  D(1,N)=X<I+1,J) 

3060  E(1,N)=R(IM,J) 

3070  F< 1 ,N) =M ( I  H , J ) 

3080  G( 1 i N) =T( 1*1 , J) 

3090  GOTO  2080 

3130  PRINT  "TRANSONIC  THRQAT  FLOW  AND  SUPERSONIC  CONICAL  NOZZLE,  1/12/ 
•3140  PRINT  "KLIE5EL  THROAT  SOL U" ION" 

3150  PRIMT  "USE  OF  KLIEGEL  S  THIRD  ORDER  THROAT  FLOW  SOLUTION  " 

5160  REM:  CONVERGENCE  EXPECTED  FOR  R=1 
5170  PRINT  "ENTER  GAMMA="; 

5180  INPUT  GO 
5190  PRINT  "GAMMA="GO 
5200  PRINT 
5210  A  ♦  =  "  N  0 " 

5220  PRINT  "ENTER  EXIT  RADIUS1"; 

5230  INPUT  00 

5240  PRINT  "EXIT  RAPIUS="0O 

5242  PRINT  "ENTER  X/R,E  FOR  FLOW  PROFILE  DETERMINATION"; 

5243  INPUT  LO 

5244  PRINT  "FLOU  PROFILE  DETERMINED  AT  >7R,E-"L0 
5250  GOTO  5290 

5260  PRINT  "ENTER  EXIT  MACH  ft  =  "; 

5270  INPUT  MO 

5280  PRINT  "EXIT  MACH  4=”M0 
52«0  PRINT 

5300  PRINT  "ENTER  THROAT  RADIUS  OF  CURVATURE*"; 

’7 3 1 0  INOIJT  CO 

5320  PR!  IT  "THROAT  RADIUS  OF  ClK'"A TORE -"CO 

5330  PRINT  "ENTER  CONICAL  WALL  DIVERGENCE  ANGLE  (PEG)5"; 

5340  INPUT  TO 

5350  PRINT  "UALl  DIVERGENCE  ANGLE="TO" ( DEG ) " 

5360  PRINT  "ENTER  ft  OF  POINTS  FOR  NET=": 

5370  INPUT  K 

5380  PRINT  "NUMBER  OF  INPUT  POINTS*";* 

5390  10=70/57.2958 
5400  X7=C0*SIN'T0) 

5410  R?= 1+CO  * ( 1  - SQR ( 1 - ( X 7/CO ) * 2 ) ) 

5420  PRINT 
5430  GOTO  5440 

5440  PRINT  "NUMBER  OF  INITIAL  POINTS="K 
5450  FOR  X=0  TO  X7  STEP  ( X7/1 0 ) -1 E-06 

5460  PRINT 

5470  Y0=1 +C0*( 1 -SQR( 1  - ( X/CO ) *2) ) 


5480  PRINT  "X="X 
5490  J=t 

5500  FOR  Y=0  TO  TO  STEP  (Y0/(K-1)) 

55! 0  OOSUB  1660 
5520  P ( 1 ,  J )  =  X 
5530  Q( t , J)=Y 
55i0  J=J+1 
5550  NE'-'T  Y 

5560  IF  B ( I , 1 ) > ! . 025  THEN  5590 

5570  NEXT  X 

5580  GOTO  5440 

5590  FOR  J=K  TO  1  STEP  (-1) 

5600  X( 1 , K+1 - J) =P( 1 , J ) 

5610  R(1,K+1-J)=G( 1, J) 

5620  M< 1 ,K+1 -J)=B< 1 , J) 

5630  T ( 1 f K+1 - J ) =C< 1 ,  J ) 

5640  NEXT  J 

5650  GOTO  3! 0 

5660  Z0=X*SQR(2*C0/(G0+-1 )) 

5670  L 1  =  Y * Y/2-0 . 25  +  20 
5680  L2=<(2*GO+9)*Y‘4-(4*GO+15)*Y*Y)/24 
5690  L2=L2+ ( 1 0*G0+ 57 ) /233+Z0+ ( Y  +  Y  -0 . 625 )  -  ( 2  *  GO  -  3 )  *  Z  0 ' 3  /  6 
5700  L3  =  ( 556  +G0"'2+1 899*00+3231  )*Y"  6/ 10368- (388*GO*2+1  233  »=G0 + 1  933 )  *  Y "  4 / 230  4 
5710  L 3 = L 3  + : 3 0 4  *  G 0 '  2  +  858  -50  + 1 2 6 9  > *Y*Y/1 728- < 2708 *G0 *2+7839+60+ 1421 1 ) /3294  a 
’20  (.4  =  /52*G0/'2aS'»*G0-3,5>  *Y  '  ><"33  4-(52+G0  >°Q*G0+303 )  *Y*Y/1 92 
5730  L4=l4+( 92*60  '2+ 180*60  +  639) /1  152 

5740  L5=((13*GO-27)/43-(5*GO-5)*Y*Y/8)*ZO*ZO+ZO*ZO*ZO*'  (4+G0*2-57 *60+27)/ 144) 
5750  L3=L3+Z0*L4+L5 

5760  B ( 1 , J )  =  1 +L1  / <  CO  +  1 )  +  (LI +L2)/(C0+1 ) ' 2+ (LI +2+L2+L3) (CO+1 )  "3 
5770  N1 =Y *Y*Y/4-Y/4+Y*Z0 

5780  N2=(8*G0+15>*Y‘5/72-<20*G0+45) *Y  "3/96+ ( 28+G0+75)  +  Y/288 
5790  N2=N2  +  Z0* (( 4*G0  +  9)  *Y* 3/1 2- ( 4*60+9 )  *Y/1 2 

5800  N3=(6836*G0  2  +  1 6695*G0  +  1  4211 ) * Y ’ 7 / 8294 4- ( 3380* GO "2  +8703*50+78^5) +  Y "5/1 .2824 
5310  N3=N3+ ( 3424*60*2+91 83*60+8964 ) *Y* 3/1 3824 
5820  N3=N3- ■ ’* O'  * Q 0 ' 2 + 1 9575 «30A20745 ' +Y  " S294 J 

5830  N4  =  ■ 556  *G0  "2+ 1 1 1 3  +  G0  +  981 ) *Y*5/1 728- ( 383*G0  "2+30 ! *G0+6«3)*Y "3/576 
5840  N4=N4+ ( 304 * GO ’2+6 45*60+5 49 ) *Y/864 

5850  N5=ZO*ZO*(( 52*G0*2+3*G0-33)*Y*3/1 92-( 52+G0*  2+27*G0 -? > * Y / 1  92; 

5860  N6=Z0  3* ( GO+1 ) *Y/4 

5870  N3=N3+Z0*N4+N5-N6 

5880  Z°=SQR( (GO+1 )/<2*'C0+1 ) ) ) 

5890  C(1 ,J)=Z9*(N1/(C0+1 ) ♦ <1 .5*N1 +N2 ) / ( CO+1 )  2+ ( 1 5 *N1 /S+2.5 » N2+N3 ) / (C0+ 1 ) *3) 

5900  RETURN 

5910  FOR  J=1  TO  K-S+1 

5920  X(1  , J >  =  X ( 3 , J) 

5930  R( 1 f  J) =R ! 3, J ) 

5940 

5950  TC1rJ)=T(3,J) 

5960  NEXT  J 
5970  1=1-2 
5980  RETURN 
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6060  B  =  <  0  0 + U  0  *  CO  S  <  TO )  -  Z  3 ) .-  'U  0 

6070  PR  1HT  "IS  C DC  (NOT  CO'  NOZZLE  UANTEI'?" ; 

■5080  INPUT  C$ 

6090  IF  C$C- “YES"  THEN  6120 
6100  B=FNB(B) 

6110  GOTO  6170 
6120  B=-FNB< B) 

6170  E=X4+UO*(SIN(TO)+SIN(B>) 

6175  L0=E+L0+00 

6180  GOTO  420 

6190  GOSUB  2570 

5200  GOSUB  560 

6210  FOR  J=1  TO  K-S 

6220  PRINT  I , J , X ( I , J),R(I, J) 

6230  NEXT  J 
6240  1=1+1 
6250  S=S+1 

6260  IF  K-SCO  THEN  6280 
6270  GOTO  6200 
6280  X 6  =  X (1  ,1) 

6290  R6*RO,f> 

6300  PRINT  "NOZZLE  EXIT  CHARACTERISTICS  NORMALIZED  AND  STORED  IN  FILE  SONIC;" 
6310  PRINT  "  N  X ( 1 , i )  R(1,N)  H  M  , N 1  T  ( 1 ,  N KRAD; " 

u31 1  FOR  N=I-1  TO  2  STEP  (-' ' 

6312  X ( N  + 1 , 1 I=X(N,1 ) 

6313  P  N  +  1 ,1  '=R'N,*.' 

6 3'  i  0  ’  M  + 1 , 1 )  =  1  t  *i , 1 5 
5315  T( N+1 , 1 ) ST (N, 1 ) 

6316  NEXT  N 

6217  GOTO  6500 
6320  FOR  N=i  TO  I 
6330  X(1 fN)=(X(N,1 >-XS >/R6 
6340  R ( 1 , N ) =R ( N, 1 ) /R6 
6350  .1(1  ,N.'=N(N,I ) 

5360  T( 1 ,N'=T(N,1 ! 

6370  PRINT  USING  6330, H, X' 1 , N) ,R( 1 ,N> ,M( 1 ,N; , T( 1 ,N) 

6380  :  St».  MI.MIII  nn.tmun  HS.NHstS 

6385  PRINT  H 1 , X( 1 , N) ,R( 1 , N ) f H • 1 , H) , T ; J ,N ) 

6390  NEXT  N 
6395  CLOSE  #1 
6490  STOP 
6500  X2=A6 
6510  R  2  -  A  7 
6520  M2=A9 
6530  T2=A3 
5540  X 1 =B6 
5550  RI=P7 
6560  Ml =B3 
6570  T1=B? 

5580  GOSUB  «30 
5630  X  ( 2 , 1 )  =  X  3 


<5440  R ( 2 , 1 )  =R3 
<5650  T( 2, 1 )  =T3 
4460  11(2,1 )=M3 
4670  GOTO  6320 

■4475  PRINT  "AT  CKOS 5  SECTION  X-P,E="'„0 
<5480  B  =  B  + 1 

•4685  B f  =  ( L 0 - X 1  ) / ( X 3 - X 1  ) 

■4490  V <  1 ,  B  *  =L0 
<4495  V( 2, B )=R1 ♦ ( R3-R1  )*B1 
<4700  V( 3, B)=H1 + ( M3-M1 )  »B1 
<5705  V(4,B)=T1+(T3-T1  >*B1 

<4710  PRINT  USING  4715,  B,L0,.tf  (2,B  •  ,V(3,B)  ,V'4,B)*5?.2?58 
■471 5  :»»»»*#».  #«.####  MH.MIHH*  RH.HHSI1  35*.H#|»S 

<4720  RETURN 

<4730  PRINT  "AT  CROSS  SECTION  X/R,E="L0 
<4735  B=B+1 

■4740  B2=(L0-X2)/ ( X 3 - X 2 ) 

<4745  V ( 1 , B )  - L 0 
<5750  V(2,B)=R2+(R3-R2)+B2 
.5755  V(3,B)=H2+(H3-M2)*B2 
4760  V(4,B)=T2+(T3-T2)*B1 

<•5765  PRINT  USING  5715, 3, L? ,  V  !  2 ,  B  ’> , 3,  B ' ,  7  ( 4 ,  B )  *57, 2'?cj8 
4770  RETURN 


Program  JETC 


Program  JETC  calculates  the  inviscid  plume  shape  using  the  MOC. 
The  last  characteristic  line,  calculated  and  stored  in  file  "SONIC"  by 
program  CDCSECT,  serves  as  a  starting  condition. 


i  nput 

quantity 

symbol 

expl anation 

y 

GO 

specific  heat  ratio 

MO 

exit  lip  Mach  number 

re/r* 

00 

exit  radius 

MJS 

M6 

jet  surface  Mach  number 

K 

number  of  points  in  MOC  net 

output 

x/r'e1',  r/re 

and  9  (deg) 

for  the  plume 

tAl 1  plume  dimensions  are  non-dimensi onal i zed  by  the  exit  radius. 


00005  PRINT"PROGRAM  JETC" 

OOOIO  PRINT  "JET  BOUNDARY  FOR  GIVEN  NOZZLE  EXIT  CONDITIONS, HHK f  1/31/78 
00015  PRINT  "FLOU  THROUGH  NON  IDEAL  C.D  NOZZLE" 

00020  DIH  A$(3),B$<3) 

00025  PRINT  "ENTER  GAHHA=" ; 

00030  INPUT  GO 
00035  PRINT  "GAMMA="GO 
00040  PRINT 

'.,0045  PRINT  "15  EXIT  MACH  H  GIVEN?": 

•  *005  ■  INPUT  At 
00055  IF  >jr’=  1  »E'j"  -0 

i.V'6’)  PRINT  •' t « T E R  E X I  T  A 0 1 1 J 5  = " ; 

00065  INPUT  00 

00070  PRINT  "EXIT  RADIUS="00 
00075  GOTO  110 

00080  PRINT  "ENTER  EXIT  MACH  #="; 

00085  INPUT  MO 

00090  PRINT  "EXIT  MACH  N="HO 

00095  PRINT  "ENTER  JET  SURFACE  MACH  M"; 

OOIOO  INPUT  M6 

00105  PRINT  "JET  SURFACE  MACH  H="H4 
00110  PRINT 

00115  PRINT  "ENTER  #  OF  POINTS  FOR  NET=" ; 

00120  INPUT  K 
00125  KO=K 

00130  PRINT  "NUMBER  OF  INPUT  F'OINTS=";K 
00135  PRINT  " INPUT  FROM  SONIC  FILE?"; 

00136  INPUT  B$ 

00137  IF  B$="YES"  THEN  1390 

00140  DIM  X ( 2 1 , 1 1 ) , R ( 2 1 , 1 1 ) , M ( 2 1 , 1 1 ) , T ( 2 1 , 1 1 ) 

00145  DEF  FNA( X ) =ATN( 1 /SQR( X*X- 1 ) ) 

00150  DEF  FNS(X)=SQR(2)*X/SGR<G"  +  ;  - 1  ■  -  <  -  X ) 

00155  iZ-  :  '*  ‘i  <  -  !  :  /;  ft  ‘  •> ' .) 

"01  60  F«FiA'  =  3^'T)tSIN(A)/SIN<TfA) 

00165  DEF  FNG(A)=SIN(T)*SIN( A)/SIN(T-A) 

001?0  PRINT  "INITIAL  POINTS  ALONG  CHARACTERISTIC  CURVE" 

00180  1=1 
00185  PRINT  "1=1" 

00190  FOR  J=1  TO  K 

00195  INPUT  X(l,J),R(1,J),M(1,J),T(1rJ) 

00200  PRINT  J,X(I,J);R(I,J);M(I,J);T(I,J) 

00205  NEXT  J 
00210  T 0  =  T (1,1) 

00215  PRINT 
00220  PRINT 

0022 5  PRINT  "FLQUFIELD  COMPUTED,  ALL  POINTS  PRINTED" 

00230  PRINT"  J  X  R  M*  THETA" 

00235  PRINT 


00240  PRINT 
00245  S=1 

00255  IF  I>1  THEN  270 
00260  GOTO  850 
00265  1=1 
00270  GOSUB  285 
00280  GOTO  365 
00295  FOR  J=l  TO  K~i 

•  )02°0  '  1  =  >  <  I  + 1  „  j  • 

•  i  ■) o "  <  2  =  i+i  \ 

'"■ICO  P1=R<I+1..” 

'0305  &2=R'I,J+1  .■ 

00310  H1=M(  IH,J) 

00315  H2«N(I,J+1) 

00320  T1 =T( 1+1 , J) 

00325  T2=T ( I, J+l ) 

00330  GOSUB  395 
00335  X ( 1+1 ,J+1 )=X3 
00340  R ( 1+ 1 , J+1 ) =R3 
00345  H( 1+1 , J+1 )=M3 
00350  T( 1+1 , J+1 ) — T3 
00355  NEXT  J 
00360  RETURN 
00365  GOSUB  635 
00370  PRINT 
00375  1=1+1 
00380  S=S+1 

00385  IF  IXN-1 )  THEN  1210 
00390  GOTO  255 
00395  H9=1 .01 
00400  A1 =FNA(FNS(M1 ) ) 

00405  A2=FNA(FNS(H2) ) 

00410  A4-A1 
00415  A5=A2 
<>0420  T4=T1 
00425  T5=T2 
00430  M 4 — HI 
00435  N5*M2 
00440  R4=RI 
00445  R5*R2 

00450  X3S  < R2-R1 +X1  *TAN(T4-A4)-X2*TAN< T5+A5 > )/( TAi1(T4  -  A4’ -T,iN< T5+m5)  ) 
00455  R3=R1 ♦ ( X3-X1 ) *TAN( T4-A4 ) 

00460  T=T4 
00465  A=A4 

00470  f*  =  r»4 

•  =PNQ(<e-  • 

•  •  .*  A  5 1  *  }  i  -  C  *j r.  ^  \ 

0 0 4 : 5  T-T5 
,;0«90  6  = 

■<0495  rt  =  l5 
>050'*  Q2=fN0<A) 


00505  IF  R2=0  THEN  665 
00510  F2=FNF(A) 

00515  H3=(m-T2)+GI*<R3-R1)/R4+F2*(R3-R2)/R5+Ql*M1+Q2*M2)./< 01+02) 

00520  T7=Q2*(M3-H2)-F2*(R3-R2)/R5 
00525  T3=T2+T7 

00530  IF  ABS<N9-H3)<0.0001  THEN  580 

00535  M9=M3 

00540  T4=(T1+T3)/2 

00545  T5=(T2+T3)/2 

00550  M4=(H1+M3)/2 

00555  N5=(H2+H3)/2 

00560  A3=FNA(FNS(M3) ) 

00565  A4=(A1+A3)/2 
00570  A5=( A2+A3)/2 
00572  R5=(R2+R3)/2 
00573  R4=(R1+R3)/2 
00575  GOTO  450 
00580  RETURN 

00585  FOR  J=K-S+1  TO  1  ‘STEP  (-1) 

00590  X ( 1+1 , J+1 )=X(I+1 , J) 

00595  R( 1+1 , J+1 )=R< 1+1 , J) 

03600  M ( 1  +  1  f  J+1  )=MCI+1 ,J> 

03605  T( 1+1 , J+1 )=T( tH , J; 

•  1 0  f- !  +F  YT  j 
006 15 
'•••>20 

•0625  IF  I-’  T  'EN  910 

•  '0633  GOTO  365 

00635  PRINT  "ROW  #=’*1+1 
00640  FOR  J=1  TO  K 

00645  PRINT  USING  650, J,X (I +1 , J) ,R< 1+1 , J ) , H (1+1 , J) , T< I +1 , J ) *1 80/3. 1 41 5? 
00650  :  J=«.  Rs*#.«ltft«  «=*#.# HR*  T=»N.NII#R 

00655  NEXT  J 
00660  RETURN 

00665  M3=< Q2+M2/2+T1 +Q1 *H1 +G1 *(R3-R1 ) /R4 )/( 02/2+81 ) 

00670  T7=G2*(H3-N2)/2 
00675  GOTO  525 
00680  XI =X ( I , K-S+1 ) 

00685  R1=R(I, K-S+1) 

00690  HI =H( I , K-S+1 ) 

00695  T1  =  T( I , K-S+1 ) 

00700  A1=FNA(FNS(M1 ) ) 

00705  H9=1 .01 
00710  A 4 = A 1 
00715  T4-T1/2 
00720  H4=M1 
00725  R4-R1/2 
00730  T9=TANX A4-T4) 

00735  X3=X1 +R1/T9 
00740  G1=FNG(A) 

00745  Q1=FN0(A) 


O0759  *3=N1+<T1*3.  mis?  '’90-0  '  .•  / Q i 
00755  IF  ABS(N3-H9H0.0001  THEN  785 

00760  N9=N3 

00765  M4=(M1+H3)/2 

00770  A3=FNA(FNS<H3)) 

00775  A4= ( A1 +A3 )/2 
00790  GOTO  730 
00785  X( 1+1  ,K-S+1 ) - X 3 
00790  R( 1+1 ,K-S+1 >=0 
00795  M ( I ♦ 1 ,K-S+1 )=M3 
00800  T( 1+1 , K-S+1 )=0 
00805  RETURN 
00810  FOR  J=1  TO  K-S 
00815  X(1 ,J)=X(3,J) 

00820  R(1,J)=R(3,J) 

00825  H(1,J)=M(3fJ) 

00830  T < 1 , J ) =T  < 3 , J) 

00835  NEXT  J 

00840  1=1-2 

00845  GOTO  630 

00850  L0=SQR( (GO+1 )/(GQ-1 ) ) 

00855  PRINT 

00860  PRINT  "PRANDTL-HEYER  EXPANSION* 

00865  PRINT 

00870  DEF  FNB(X)=X*SQR(  (GO+1  )/(2+(G0- !  >*X*X) ) 

03875  DEF  FNC( X )=SQR( (X*X-1 )/( (GO+1 ) / (GO- I ! -X*X>) 

00880  DEF  FN0(X)=L0*ATN(F.v<:-  *'  -  ’-AT* l L$*F»r  r, ,  \ 

•  n-  995  i  ■-  =  (7 m.i  Priii -6 )  >-FN0  >  F NB  '■  HO '  >  +  TO  .  * 1 80/3.1  4 1  59 

•0990  PRINT 

00895  PRINT  "INITIAL  SLOPE  ( DEG ) =ML9 

00900  1=1 

00905  P=2 

00910  H ( 1 )  =0 

00915  I ( 1 )  =  1 

00920  J(1)«L9*3. 14159/130 
00925  J=1 

00930  PRINT  "ENTER  •  STEF‘S="; 

00935  INPUT  N 
00940  N0=N 

00945  PRINT  "»  OF  5TEPS=”N 

00950  FOR  H=HO  TO  H6  STEP  (M6-M0)/(N-J ) 

00955  Ml =FNB(H) 

00960  0=FN0(H! ) 

00965  PRINT  "H=KH; "M*=*M1 ;*OHEGA="0; "DEG="0*1 80/3. 1 41 5? 
00970  X( I , J) =0 
00975  R< I , J)=1 
00980  M(I,J)=H1 

00985  T(I,J)=0-FNO(FNB(H0))+T0 
00990  1*1+1 
00995  NEXT  H 
OIOOO  GOTO  265 


01005  J=1 
01010  X2=X(N, J+1 ) 

01015  R2SR(N, J+1 ) 

01020  M2=M < Nf J+1 ) 

01025  T2=T(N,  J+'l ) 

01030  A2=FNA(f'NS 1 *  l  ■  > 

01035  A3=PNA'FNS(INN,  J  1  1  * 

01040  R!=R':N,J> 

01045  X1=X(N,J> 

01050  Ml =H ( N , J ) 

01055  T1=T (N,J) 

01060  A5=A2 
01065  M3=FNB(H6) 

01070  A3=FNA(FNS(M3) ) 

01075  T5=T2 
01080  M5=M2 
01085  R5=R2 
01090  T9=T1 
01095  X3=X2 

01 100  X3=(R1 -R2+X2+TAN(T3+A5)-X1 *TAN(T9) )/(TAN(T5+A5) -TAN(T9) ) 

01105  R3=R2+(X3-X2)*TAN( T5+A5) 

OHIO  Q2- 1  /  ( M5*T  AN  <  A5 ) ) 

01115  F2=SIN(A5)*SIN(T5)/SIN(A5+T5) 

01120  T3=T2+Q2*(H1 -N2)-F2*(R3-R2)/R5 
01125  T9  =  < T 1 +T3 ) /2 

01130  X7=( R1 -R2+X2*TAN( T5+A5) -XI *T AN < T9 ) )/(TAN(T5+A5)-TAN(T9) ) 

01135  IF  ABS(X3-X7)<0.001  THEN  1170 

01140  X3=X7 

01145  T5= ( T2+T3 ) /2 

01150  M5=(M2+M3) /2 

01155  A5=(A2+A3)/ 2 

f>11 6  0  95  =  f  R2+R3  '  "2 

•ni65  GOTO  •  105 

01170  RO  = '  !*  1 -R3 ) .  ' C C3 <  T 1 )  -COS  ( T3 ) ) 

01! 75  PRINT  "RO="RO; "X3="X7; "R3="R3; "T3 (DEG )="T3+1 80/3 . 1 4159 J "M3-"M6 

01180  X<N+1 ,1 >=H(P)=X3 

01185  K(P)=R0 

01190  R( N+1 , 1 ) =1 ( P ) =R3 

01195  H(N+1 , 1 )=M3 

01200  T (N+1 , 1 )=J(P1=T3 

01205  RETURN 

01210  J=1 

01215  GOSUB  1005 

01220  GOSUB  1250 

01225  N=N+1 

01230  P*P+1 

01235  IF  N'NO+KO-2  THEN  1330 
01240  K=K-1 
01245  GOTO  1215 
01250  FOR  J=1  TO  K-2 
01255  X 1 =  X ( N + 1  , J) 


0)260  X2=X(N, J+2) 

0)265  R1=R(N+1,J) 

0)270  R2-R < N , J+2 ) 

01275  M1=H(N+1,J) 

01280  M2=H<N, J+2) 

01285  T 1 =T ( N  +  1 , J ) 

01290  T2=T(N,JO 
01295  G0SU8  395 
0)300  X(N+1 ,J+1 )=X3 
0)305  R ( N+1 , J  +  1 ) =R3 
0)310  H ■  • 1  ♦ * ,  > i  *s03 
01: '5  '•*»*),  :-!  '*i"» 

■‘320  NEYT  J 
)1325  F' E TIJ P N 

01330  PRINT  "''•TER  5EGNENT  #*"; 

01335  INPUT  P 

01340  P«P+1 

01345  PRINT  "P="P 

01350  PRINT  "INPUT  RADIUS1"; 

01355  INPUT  Q 
01360  U=J(P-I) 

01365  U)=C0S1U) 

01370  U2=U1 +(Q-I(P~1 ))/K(P) 

01375  U 2 = A T N  <  S Q R  <  1  - ;J 2  U 2 )  / U 2 > 

01380  PRINT  "AT  RADIUS=,’Q"THETA=,,U2 +  180 /3. 1 41 59 ; “X*"H(P- ! ) +K ( P ) *( SIN (U ) -SIN 1  U2 ' ) 
01385  GOTO  t  330 
01390  1=1 

01400  FILE  HI ="SGNIC" 

01405  FOR  N=)  TO  K 

01410  INPUT  »1  ,D<1  ,N),E<1,N>,F'1,N),GU,N) 

01415  NEXT  N 
01420  FOR  J=1  TO  K 
01425  X(1,J)*D(1,J) 

01430  R(1, J>*E<1,j) 

01435  H(‘,J'SFM,J) 

01440  TM  ,J.)  ,  J) 

01445  NEXT  J 
01450  got:  ; 


57 


Program  NOCON 

Program  NOCON  uses  a  second-order,  non-ideal  conical  source  flow 
approximation  to  determine  a  plume  shape.  The  reason  for  running  this 
program,  besides  as  a  comparison  to  Program  JETC,  is  that  as  an  output 
the  initial  deflection  angle,  9p,  and  the  initial  radius  of  curvature, 
Rp,  are  given  (convenience).  These  parameters  serve  as  input  to  Program 
PSIM. 


input 


quantity 


symbol 


explanation 


specific  heat  ratio 
exit  lip  Mach  number 
nozzle  exit  angle  (deg) 
jet  surface  Mach  number 
measure  of  the  local  rate  of 
acceleration  near  the  nozzle 
lip  [12] 


output 


9p»  Rp»  x/rg,  r/r  ,  and  9  (deg)  for  the  plume 


250  GOTO  260 

255  C0*3*SIM<C)*SQR(FNC(B))/(B*B) 

260  Ct=C0*C2 

265  E2=Et+F*( 02-01 -FNA( 1/B)+FNA( 1 /B) > 

270  B1=(E2-E1)/10 
275  SI =S2=S3=S4=0 
280  FOR  N=1  TO  10  STEP  1 
285  E3*E1+D1*fM-1/2) 

290  SI  =C0S(E3/F> I N ( E 3 )  * (-0.5) *C0S(E3) *  ( -1  /( 2 :*=  F  '' 2 ) )  *:B1  *31 
295  S2sSIN(E3/F) +SIN(E3) * ( -0.5)*C0S(E3) ' ( -1 / (2*F‘2) )*EH  +S2 
500  S3=SIN ( E3/F )  *SIN( E3 )  ( - 1 .5) *COS  (E3K  ( <  A-3)/ (2*  <  A-1  ) )  > *D1  +S.3 
.505  S4=C0S(E3/F)*SIN(E3)  (-1 .5)  tC0S(E3)'  (  !A-3)/'2MA-1 )  >  »D  i+24 
"5 1 0  NEXT  N 

315  U0  =  -  ( ( COS  ( E2  >  "  ( ( 3  t:A  - 1  ?  /  ( 2*<  A-1  1 ) ) ) * S I N «' E 2  ■  "0.5)  2‘F' 

•  20  «'  (31  -ir*S3'-  *COS(BO>  +  ' S2  +  -*34»  *SIN(B0)+C1 ) 

.525  C2=FNA!  1/D) 

530  T2=C+Q2-01 

,335  R0=-  ( ( U1 /SQR ( J) +SIN( T2)*SIN(G2 > *2) / (SIN* 2*G2> ) ) 

.540  R 1  =  - 1  /  ?:  0 
.345  PRINT 

350  PRINT  "SOLUTION" 

.352  T2  =  T2*57.296 

.355  PRINT  "INITIAL  SLOPE  OF  PLUME  THETA-F="’2 
360  PRINT  "(INITIAL)  RADIUS  OF  CURVATURE  R 1  ~ "  R 1 
365  PRINT 
1370  PRINT 

375  PRINT  "PLUME  SHAPE  APPROXIMATED  BY  CIRCULAR  ARC" 

380  T2=T2/57. 296 

390  FOR  S=0  TO  1.5  STEP  0.1 

395  T  =  FNA(SIN(T2)-(S/R1) ) 

400  S=R1*(SIN(T2)-SIN(T) ) 

405  R=J+RJ»(C03(T? -C0S(T2) ) 

110  PRINT  s ; R;  T *57.2'56 

'20  NEXT  S 

125  PRINT 

430  PRINT 

425  END 

READY. 


Program  WSHOCK 


Program  WSHOCK  uses  Weak  Shock  Modeling  to  calculate  an  equivalent 
jet  surface  Mach  number  to  model  a  prototype  gas  with  a  different  spe¬ 
cific  heat  ratio  having  the  same  pressure  rise-deflection  angle  re¬ 
lationship. 


input 

quantity 

YP 


symbol 

K1 


Ml 


K2 


explanation 

specific  heat  ratio  of  the  prototype 
gas 

jet  surface  Mach  number  for  the 
prototype  gas 

specific  heat  ratio  of  the  model 
gas 


output 

Mpm  M2  jet  surface  Mach  number  for  the 

model  gas 


5  PRINT  "WEAK  SHOCK  MODELING,  JET  SURFACE  MACH  NUMBERS,  HHK,1 ,15,80" 

10  PRINT  "ENTER  GANNA, PROT.  K1  AND  GANNA  MODEL,  K2="; 

20  INPUT  K 1 , K 2  3 

30  PRINT  "GANNA,  PR0T0TYPE="K1 ; "GANNA ,  H0DEL="K2 
40  PRINT  "ENTER  N,F  JET  MACH  #,  PROTOTYPE2"; 

50  INPUT  N 1 
60  PRINT 

70  A=M1 *4/(N1*N1 -I ) 

(30  B=A*K1*K1./(2*K2*K2) 

90  X=B*SQR(B*B-A+K1*K1/(K2*K2>) 

100  N2=SQR(X) 

110  PRINT  "M,F-PR0TOTYPE="M1 ; "N , F -MODEL  ="M2 
120  GOTO  40 

130  END  £ 

READY . 


Program  SECTP 


Program  SFCTP  calculates  flow  properties  at  a  given  cross  section 


within  the  plume  using  the  MOC. 


input 


quantity  symbol 


expl anation 


specific  heat  ratio 


jet  surface  Mach  number 


number  of  points  in  the  MOC 


cross  section  at  which  flow 


properties  are  to  be  calculated 


x/r  ,  r/r  ,  M*  and  9  (deg)  characteri sties  stored  in  file 


'SONIC' 


output 


x/r  ,  r/r  M*  and  9  (deg)  at  the  prescribed  cross  section 


1 


M 


\  .** 

y:-: 

r 

ErV. 


IS 


64 


5  p; 3  i"  r»" 

!5  PRINT  "rLOU  T»':;'j;H  C-'C  MZI'.S,  ErrrECT  OF  ”^'C-T  :.  ’  .  f ' ,  :..■.»  •  =  - 1  -*fts  ■  r 

ift  PRINT  " PL'J4E  -L'C’J  CALCULATED  FQF:  GIVEN  ACROSS  SECT!  j ■-  ' 

22  FILE  81 = "SONIC" 

25  PIH  A«(3),B«(3),C*(3) 

27  Q0=3. 1 41592 

30  PRINT  "ENTER  GAMMA="; 

35  INPUT  GO 

40  PRINT  "GAMMA="GO 

45  DEF  FNE(X)»ATH(X/SSR(1  -X*y«-fE-??)  * 

30  PRINT 

105  PRINT  "ENTER  JET  SURFACE  MACH  IN"; 

NO  INPUT  H6 

115  PRINT  "JET  SURFACE  MACH  »=“M6 
120  Y=0 
125  PRINT 

130  PRINT  "ENTER  8  OF  POINTS  FOR  NET-" ; 

135  INPUT  K 
140  K 0 = K 

145  PRINT  "NUMBER  OF  INPUT  POINTS*" ;K 

155  DIM  Xmf11),R(11,11),N(11,11),T<11f11),A<4,2l)fV(4f11> 

160  PRINT  "ENTER  Z0>0  IF  FLOU  CROSS  SECTION  TO  BE  CALCULATED  AT  ZO"; 

165  INPUT  ZO 

170  PRINT  "FLOU  CROSS  SECTION  DETERMINED  FOR  ZO="ZO 
175  DEF  FNA (X)=ATN( 1 /SQR(X*X-1 ) ) 

180  DEF  FNS ( X )  =  SQR  T 2 ) * X / SQ R < GO  + 1  - ( G 0  - 1 ) *X*X ) 

135  DEF  FNQ ( A ) 8 1 / C  M *  TAN ( A ) ) 

190  DEF  FNF(A)=SIN(T)tSIN- A) -Ilri'TiA) 

195  DEF  FNG ( A ) =S IN ( T 1 *S I N ( A ) / SIN ( T-A ) 

200  PRINT  "INITIAL  POINTS  ALONG  CHARACTERISTIC  CURVE" 

210  1=1 

215  PRINT  "1=1" 

220  PRINT  "DATA  FROM  FILE=Of FROM  KEY*1"; 

225  INPUT  L 

230  IF  L=0  THEN  1 365 

235  FOR  J=1  TO  K 

237  PRINT"ENTER  X/RE,  R.'RE,  M«,  THETA (RAD)  FOR  POINT  I  "J 
240  INPUT  X(1, J> ,F( 1 , J),M(1 ,J)fT(t,J> 

245  PRINTJ,X(I,J),R(I,J},N(If J>,T(I, j) 

250  A< 1 f J ) =X ( 1 , J) 

255  A(2,J)=R(1,J) 

260  A ( 3, J) =M( 1 , J) 

265  A(4,J)=T( I ,J) 

270  NEXT  J 

275  PRINT  "SHALL  DATA  IE  STORED  IN  FILE  1 (SONIC)"; 

280  INPUT  B* 

285  IF  B$<>"YES"  THEN  295 
288  FOR  J=1  TO  K 


v 


3 


■.  •% 


A 


.  t 

Zt 

J 


V 


8 

<T 


y 

*  d 


y 


V 


J 


m  - 

I 


* 


5 


. •  \v, 


28 9  PRINT  #1,  X(?,J)IR(1,  J>,nn, J\TM,  J) 

290  MEXT  J 
292  CLOSE  #1 
295  TO = T (1,1) 

300  PRINT 

102  M0=FNS(A<3, 1 ) ) 

303  PRINT"EXIT  MACH  S  =  "HO 
305  PRINT 

310  PRINT  "FLOUFIELD  CONFUTED,  ALL  POINTS  PRINTED 
315  PRINT"  J  X/RE=  R7RE 

320  PRINT 
325  PRINT 
330  S=1 

340  IF  IM  THEN  355 

345  GOTO  965 

350  1=1 

355  GOSUB  365 

360  GOTO  475 

365  FOR  J=l  TO  K-1 

370  X1=X( 1+1 , J) 

375  X2=X ( I , J+ 1 ) 

380  R1 =R« 1+ I , J) 

385  R2=R( 

390  HI aH( I *1 , J) 

395  H2=H(I,J*1) 

100  T1«T(IH,J) 

105  T2»T(I,JM) 

110  GOSUB  535 

115  Xa  +  1,JM)=X3 

*120  R(I+1 , J+1  )=F'3 

125  H( 1+1 , J+I )=M3 

130  T'  I*t ,  J«-i  )~-*3 

135  Ir  X3>"0  AND  V  -20  THE'-'  145 

140  GOTO  450 

145  GOSUB  1640 

150  IF  X3: 20  AND  X2CZ0  THEN  460 

155  GOTO  465 

460  GOSUB  1635 

465  NEXT  J 

170  RETURN 

175  GOSUB  745 

180  PRINT 

485  I  =  IM 

490  S=S+1 

495  IF  IXN-1)  THEN  1360 
500  GOTO  340 
505  N9= 1.01 
510  A1=FNA(FNS(H1 ) ) 

515  A2=FNA(FNS(M2) ! 

520  A4=A1 


530  T 4*  T I 
535  T5*T2 
540  N4=M1 
545  H5=M2 
550  R4=R1 
555  R5=R2 

560  X3*(R2-R1+X1*TAN(T4-A4)-X2+TAN<T5+A5>)/(TAN(T4-A1)-TAN<T5+A5>) 

545  R3=R1+(X3-X1 >*TAN(T4-A4) 

570  T=T4 
575  A=A4 
580  N=H4 
585  Q1=FNQ(A) 

590  G1 =FNG ( A ) 

595  T=T5 
600  A=A5 
605  M=«5 
610  Q2=FNQ'A> 

6 1 5  IF  R2=0  THEN  ’SO 
620  F2=FNF ( A ' 

625  M3= ( (T1 -T2) *80/ 1 90+61  * < E3-R’ • /R4  +f 2* < 93-R2 3 /P5+91 »M1 +82*12)/ ( Q 1 +02) 
630  T7«Q2*(N3-H2)-F2*(R3-R2)/R5 
635  T3=T2+T7 

640  IF  ABS(H9-H3)  '0.0001  THEN  690 

645  H9  =  M3 

650  T4=(T1 +T3)/2 

655  T5=(T2+T3)/2 

660  M4=(M1  +  M 3 ) / 2 

665  «5=(H2+M3)/2 

670  A3=FNA(FNS(M3) ) 

675  A4=(A1+A3)/2 
680  A5=( A2+A3)/2 
685  GOTO  560 
690  RETURN 

695  FOR  J=K-S+1  TO  1  STEP  <-1) 

700  X(I+1, J+1)=X(I+t,J) 

705  R ( I + 1 ,  J+1 )=R(I+1 ,J) 

710  H(I+1,J+1)=M(I+1,J) 

715  T< I  +  T , J  +  1 >  =T( 1  +  1 , J) 

720  NEXT  J 

725  K=K+2 

730  GOTO  715 

735  IF  1=2  THEN  925 

740  GOTO  475 

745  PRINT  "RIGHT  RUNNING  CHARACTERISTIC  »"I+I 
755  FOR  J=1  TO  K 

760  PRINT  J,X(  1  +  1,  J),R(I  +  i,J),f1'I  +  »,  J),T<I  +  1,J>  *18O/Q0 
7 70  NEXT  J 
775  RETURN 

780  H3=  (Q2*M2/2  +  T1  *00/ 1 80  +  01  +M1+G1  +={R3  -R1 ) /R4 ) / ( Q2/2+Q 1  ) 

785  T7=Q2*(M3-M2)/2 
790  GOTO  635 


X 1 =  X ( I ,  K-S+1 ) 

R1=R(I, K-S+1 ) 

M1=h(I, K-S+1) 

T1=T(I, K-S+1) 

At  =FNA(FNS<M1 )) 

M9=1 .01 

A4=A1 

T4=T1/2 

«4=M1 

R4=Rt/2 

T9=TAN(A4-T4) 

X3=X1+R1/T9 
G1 =FNG(A4) 

Q1=FNQ(A4) 

N3=H1+(T1 *00/180-61 )/Q1 

IF  ABS(M3-M9K0.0001  THEM  900 

M9=M3 

M4=(N1 +M3)/2 
A3=FNA(FNS<H3> ) 

A4=(A1 +  A  3  > / 2 

GOTO  345 

X'I  +  1 ,K  3+1 ) =X2 

R(I+1,K-S+1)-0 

MC  +  t  , K-S+1  )-M3 

T(  1  +  1 , K-S+1 )  =0 

RETURN 

fq '>  jii  rn  «;-? 

X!', J'=X'3,J> 
pm,j'=r:3,J) 

T(  1 ,  J)  =  T<3, J) 

NEXT  J 
1*1-2 
GOTO  740 

L0=3GR< (GO+f )/<30-1 > ) 

PRINT 

PRINT  "PRANBTL-MEYER  EXPANSION" 

PRINT 

DEF  FNB(X)=X*SQR( (GO+1 )/(2+ (GO- I) +X *X ) ) 
BEF  FNC(X)=SQR(  (X*X-1 )/( (GO+1  )/(GO-l  ?-X*.X 
DEF  FNO(X)=LO*ATN(FNC(X))-ATN< LO+FNC (X ) ) 

)  L9= (FNO(FNB(M£) ) -FNO(FNB(HO ) ) +T0) *  I 60/00 
5  PRINT 

)  PRINT  "INITIAL  SLOPE  (BEG)="L9 
5  1  =  1 
)  F  =  2 
5  H ( 1 ) =0 
J  I(1)=1 

5  J(1)  =1.9*00./!  30 
}  J =1 


INPUT  N 
NO=N 

PRINT  *N  OF  STEPS="N 

FOR  M=MO  TO  Ho+1E-05  STEP  (M6-M0 '/<N-1 > 

Ml =FNB!M) 

0=FN0(M1 ) 

PRINT  ,,M=,,M;,,M*  =  ,,M1 ; "QMEGA="0 ; "t  £3="0+1 3O/Q0 
X ( I , J ) =0 
R(I,J)=1 
M(I, J)=M1 

T( I , J) =0-FN0( FNB (NO ) ) +T0 
1  =  1  +  1 
NEXT  M 
GOTO  350 
J=1 

X2=X(N, J+1 ) 

R2=R(N, J+1 ) 

M2=M(N, J+1 ) 

T2=T (N, J+1 ) 

A2=FNA(FN3(M2) ) 

A3=FNA(FNS(M(N, J) )) 

R1=R(N,J) 

XI  =X  <N,  J ) 

M 1  =  M  ( N ,  J ) 

T1 =T (N, J) 

A5=A2 

H3=FNB(«6) 

A3=FNA(FNS(M3> ) 

T5=T2 
M5-H2 
R  3=R2 
T<?=T1 

xj=x: 

X3=  <R1 -R2+X2*TAN(T5 :A5 ) -X 1  +TAN( T?) ) /  ( TAN  •  TT+A5 .'-"AMI T?  ■ ' 
R2  =  P2  +  ( X 7 - X 2  kTAN^T=j-45' 

Q2  =  1/M5*TAN(~5> ) 

F2=SIN( A5) +3IN(T5)/SIN(A5+T5> 
T3*T2+Q2*(H1-M2)-F2*(R3-R2)/RS 
T  9  = ( T 1 +T3)/2 

X7« ( R1 -R2+X2  *TA»  <  T5+A5 ) -XI *TAN  <  T?'  > .’  < TAN<  T5+A5 ) -TAN ( T? ) » 

IF  ABS(X3-X?)<0.001  THEN  1285 

X3=X7 

T5  =  ( T2+T3 ) /2 
M5= ( M2+M3 ) /2 
A5=(A2+A3)/2 
R5= ( R2+R3 ) !2 
GOTO  1220 

RO=(RJ-R3)/(COS(T1 )-C0S(T3>) 

PRINT 

PRINT  "RIGHT  RUNNING  CHARACTERISTIC  #"N+ ! 

PRINT  "RO="RO 


1305  F'RINTI ,X3,R3, M3, T3  *1 90/Q0 

1310  IF  X3)Z0  AND  X2TZ0  THEN  1320 

1315  GOTO  1325 

'320  GOSUB  1730 

'  325  X(N  + V  ' -H 3 

1330 

!  335  C  f!*  1  1  '  ;  T  ■’  p  rf  ! 

:  it*  >*.t fj4.i  t)  rf^7 

13*5  '(N  +  1  ,  1  )  -=j  s  '.=’3 

1350  IF  H ( P '  Z 0  AND  H-P-I-'ZO  "!JE 

1355  RETURN 

■360  J=! 

1165  GOSUB  1 1 20 
1 370  GOSUB  1400 
!  ’75  N-'-'  +  l 
! 330  P = P *  1 

I 3S5  [F  N  NO+KO  2  THEN  1515 

1390  K*K-1 

1355  GOTO  1365 

1400  FOR  J--1  TO  K -2 

1405  X 1  =  X '  N  ♦  1 , J ) 

1410  X2=X ( N , J+2) 

■14'  5  R 1  2  R  ( N  + 1 ,  J ) 

1420  R2  =  R ( N,  J+2 ) 

1425  Ml  =M  ( N+t  ,  J ) 

1430  M2SM(N, J+2) 

1435  T1-T(N+1 , J) 

1440  T2=T(N,J+2) 

1445  GOSUB  505 

1450  X(N+1 , J+1 )»X2 

1455  R ( N+1 , J+1  )=F:3 

1440  M( N+1, J+1 ) -M3 

1465  TfN+1 , J+1 )=T3 

1  470  PRINT  j  +  1 ,  X  3 ,  R3,M3,  T3H  30/Q0 

’475  IF  X !'•••:•?  AND  XI  10  THEN  1485 

'  480  G0T0  1 490 

1485  GOSUB  ‘6*C 

:  490  If  :o  ...jfi  VG  XGO  ’HEN  1  TOO 

1465  5 1 T  o  1  e.  0  r 

1300  GOSUB  1730 

1505  NEXT  J 

1510  RETURN 

1515  PRINT  ’ENTER  SEGMENT  »■="; 

1520  INPUT  e 

1525  IF  P=0  THEN  1850 

1530  P=?+l 

1535  PRINT  "P=”F'-1 

1540  PRINT  "INPUT  RADIUS'"; 

1545  INPUT  0 
1550  U=J(P-1’< 

1555  U1=CGS(U) 
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1560 
1565 
1570 
1575 
1585 
1590 
1595 
1600 
1605 
1610 
1615 
1620 
1625 
1630 
1635 
1640 
1645 
1650 
1655 
1660 
1665 
1670 
1675 
1630 
1685 
1  690 
!  695 
1700 
1  735 
1710 
1715 
I  720 
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U2=U1 +(G-I (P-1 ) )/K(P) 

U2=ATN(SQR(1  -U2*U2)/t!2) 

PRINT  "AT  RADI US="Q"THETA="U2:*180/QO;"X="H( P-1 )+K<P)*(SIN(U) -SINUJ2) > 
GOTO  1515 
FOR  J=1  TO  K 
X( 1 , J)=A( 1 , J) 

R(1,J)*A(2,J) 

M( 1 , J)=A(3, J) 

T(1, J)=A(4,J) 

PRINTJ,X(1,J>,R(1, J>,fi(1,J),T(1,J) 

NEXT  J 

FOR  J=2  TO  K 

IF  A(?  J)  10  'i I , J  -  i ■  ’  70  THEN  ‘ 775 
NEXT  / 

GOTO  295 

PRINT  "AT  CROSS  SECTION  X-“ZO 

y=y+i 

Vt=(Z0-X1 )/( X3-X1 ) 

y  ( i ,  r )  =  z  o 

V(2,Y>=R1+(R3~R1 ) *01 
'.M3,  Y!  =H!  MM3-N1  •  +  V 1 
y  <  4 ,  Y  >  =  <  T 3 - T 1 >*V1 +T1 

P  R I N  T  Y ,  V  ‘  1 ,  Y  >  f  V (  2 ,  T ) ,  V  ( 3 ,  Y > ,  '9  ( 4 ,  Y )  M  8  0  /  Q  0 
RETURN 

PRINT  "AT  CROSS  SECTION  X="Z9 

Y  =  Y  +  1 

V2=(Z0-X2)/(X3-X2) 

V  ( 1 ,  Y )  =  Z  C 

0(2,  Y !  --R2-MR3-R2)  *  v2 
V(3fY)»N2+(H3-R2)*V2 
V(4,Y)»(T3-T2)*V2*T2 


PRINTY, 0 ( 1 , Y ) , V (2, Y ' , 9 ( 3 ,  Y ) ,' 

RETURN 


g 


4 


& 


>»l 

1730 

PRINT  "AT  CROSS  SECTION  X*":3 

£•% 

'  T35 

Y  =  Y  + 1 

4 

1740 

02-'  Z "  -  *  I .  ’  X3  -  ■•;  » 

* 

m- 

1  -jc 

y(i,  y ’ =zo 

VQ 

■i  750 

’it 2,  Y }  -PI  *  •  R3-R2:  *''2 

V* 

• 

1  755 

y(3,Y?=M2  +  (.'l3-f12)H'2 

"y~y 

1760 

V(4,Y>»(T3-T2>«V2+T2 

1 765 

PR INT  Y , V ( 1 , Y > , V ( 2 , Y » , 0 ' 3 , Y  ? , V ( 4 , Y ) « 1 80/00 

•V 

i  770 

RETURN 

1775 

PRINT  "A!  CROSS  SECTION  X  =  "ZO 

\\ 

-• 

1780 

1785 

Y  =  Y  ♦  1 

V3* ( ZO-A ( 1 f  J-1 5 )/ 1  All f  J) -A < 1 , J -1 ) ) 

V* 

(hi 

\7. 

1790 

V(1  ,Y)=ZO 

1795 

V(2,Y)=A(2,J-1  )  +  (A(2,J)-A(2,  J-1  >HV3 

.*? 

■>* 
.S  • 

1800 

V(3,Y)=A(3, J-1 )*(A(3, J)-A(3, J-1 ) )*V3 

^  * 

■S' 

,sC\ 

1 805 

V ( 4, Y )  =A( 4 , J-1 ) ♦ ( A( 4 , J 1 -A ( 4 , J  -1 ) )+V3 

1810 

PRINTY,  V(  1 ,  Y)  ,y (2,  Y)  ,0 (3,  Y)  ,V(  ■», Y) 

N.  1 

1815 

GOTO  1630 

J 

■\V 

> 

•\V. 

& 

Lt2 

J 

■g 

1820  Y=Y  +  1 

1825  V(4,Y ) =FNE <  (H(P)+K(P)*SIH(J(P))-20?/K(P)) 

1830  V(2,Y)=I(P)+J«(P)*(C0S1V(4,Y) )-C0S<J(P)> » 

1835  PRINT  "CROSS  SECTION  AT  *  =  "30 
1 840  PRI NT Y ,  30  f «  2 ,  Y )  f  N3 ,  V  '  4 ,  Y ) *  1 80/30 
1845  GOTO  1355 

1  850  PRINT  "SHALL  NEU  PLUME  BE  CALCULA TED?  =  “ ; 

1855  INPUT  Ct 
I860  IF  C$=" YES"  THEN  105 
1361  IF  C$C  "YES"  THEN  1750 
1865  PRINT 

1870  PRINT"  X/RE,  R/RE,  N:*  AND  THETA (RADI  ARE  BEING  READ  FROM  'SONIC"" 
1875  FOR  J=1  TO  K 

1  880  INPUT  #1rD'1,  J),E(1,  J)  ,F(  1 ,  J.\G(  1 ,  J) 

1885  NEXT  J 

1890  FOR  J=1  TO  K 

1  895  A(1fJ)=Xnf J)  =  D(1fJ) 

•1900  A(2,  J)=R<  1 ,  J)=E(1 ,  J) 

1905  AU,J>=NU, J)-F( l,J) 

1910  A  <  4  f  J ) =T  ( 1 ,  J  0  =G  ‘1 1 ,  J  > 

:  91 1  PRINT  A >  1 ,  J ,  A ( 2 ,  J ) ,  A •  3 ,  J )  r  A ( 4 ,  J ! 

1915  NEXT  J 
1920  GOTO  295 
1  550  END 
READY. 
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Program  PS IM 

Program  PSIM  calculates  what  flow  parameters  are  required  to 
produce  an  equivalent  plume  shape  when  a  different  specific  heat  ratio 
is  used. 


i  nput 


symbol 


explanation 


initial  plume  deflection  angle 
(deg) 

initial  plume  radius  of  curvature 
specific  heat  ratio  of  the  model  gas 
jet  surface  Mach  number  for  the 
model  gas 

trial  value  of  exit  lip  Mach  number 
for  the  model  gas 


output 

MLm  and  9^-final  values  for  the  model  yas  at  the  exit 
P-final  pressure  ratio, 

Pb/Pg,  x/re,  r/re  and  9  (deg)  of  the  model  plume 


00005  PR  IT  "  FRllGRAH  p S  IV V  [St?  :  2  ' 
f) 0 •} 20  D E r  - N A '  X )  =  A t fj { X / : 5 *  ’  - x  *<  +  :Z  5 > ) 

00025  BEF  FNB'X»=ATN<SQR(1-X*X).,X> 

0002?  DEF  FNR(N)=((2+<A-1)*H*H)/(A+t))  ( (A+1 )/<2*<A-1 ) )  )/M 
00035  PRINT  "PLUHE  SHAPE  GIVEN  ?YESO,NO  1"; 

00040  INPUT  K1 
00045  IF  K1 =0  THEN  710 
00050  Z=0 

00055  PRINT  "PROTOTYPE  FLOW  SPECIFIED" 

00060  PRINT  "ENTER  GANHA-P"; 

00065  INPUT  A 

00070  PRINT  ”GAHHA="A 

00075  PRINT  "ENTER  NOZZLE  EXIT  MACH"; 

00080  INPUT  B 

00085  PRINT  "NL-P="B 

00090  PRINT  "ENTER  NOZZLE  EXIT  ANGLE"; 

00095  INPUT  C 

OOIOO  PRINT  "THETA-L-P="C 

00102  C=C*3. 14159/190 

00105  PRINT  "ENTER  JET  SURFACE  1ACH"; 

.10110  INPUT  d 

.1 01 1 5  PRINT  "NPP  =  "0 

00120  E=? 

00135  PRINT 

00140  F-((A-1 )/(A+1))-0.5 
00145  G=FNA( 1 / B ) 

00150  H»C-G*3. 14159/2 
00155  El =ATN(F/TAN(G) ) 

00160  B0«(E1)/F-H 

001  65  BEF  FNC(X)-<  ( ( A+1  )/2)  *:X"'2)/<  1  +  ( ( A-1  )*X*2)/2) 

00170  I=FNC(B) 

00175  J=FNC(D) 

00180  X=I 
00185  GOSUB  215 
00190  01 =L 
00195  X=J 
00200  GOSUB  215 
00205  02=L 
00210  GOTO  230 

00215  K  =  < ( 1 -X)/<X-1 /F  2) )  '0.5 
00220  L=ATN<K)/F-ATN(K/F) 

00225  RETURN 

00230  C0=3*SIN<C)/(F+B) 

00235  C2=-2*F*SQR< I )/(COS (El ) " ( < ( 3  *  A ) - 1 ) / (2 M A-1 ) ) ) *SIN(E1 ) * . 5 ) 
00240  Cl  =C0*C2 

00245  E2=E1 *F* ( 02-01 -FNA ( 1 /D ) +FNA (1/B? ) 

00250  B1=(E2-E1 ) / 1 0 
00255  SI =S2=S3=S4=0 


00260  F0*  N= ‘  TO  to  STEP  t 
00265  E3=E1  +D1  KN-1/25 

■  0270  SI =C0S(E3/F)*3IN (E3) ' < -0.5)*COS (E3) " ( -1/(2*F'2) )+B1  <-'51 
00275  S2=SIN(E3/F)*SIN<E3r(-0.5)*C0S(E3>  ‘<-i/<2*F  2) )«D1+S2 
00280  S3=SIN(E3/F)*SIM(E3)'(-1 . 5 ) *COS I E3 ) * ( < A- 3 ) / < 2* ( A- 1 ) ) )*D1 +S3 
00285  S4=C0S(E3/F)*SIN(E3)'(-1 .5)  *C0S(E3K  <  <  A-3)/<2*<  A- 1 ))>*D1+S4 
00290  NEXT  N 

00295  U0=-(  (C0S(E2) ' ( (3*A- 1 )/ (2* ( A-1 ) ) ) )*SIN( E2)  0. 5)/ (2 *F) 

00300  U1=U0*( (SI -F*S3) *C0S(B0)+(S2+F*S4) *SIN(B0)+C1 > 

00305  G2=FNAC1/B) 

003’ 0  T2=C+O2-01 

0031  5  R0=- ( ( Ul/SaS  ( J )  +  5 1 M  ’  T2  » *SIN (  .52  >*2)/ !SIN(2*G2)  > ) 

•00320  Rf  =  - 1  /RO 
00325  PRINT 
■0330  ""GLUT  [ON" 

•0033 j  C'9IN T  "INITIAL  SLOPE  OF  *Ujme  THETA  '-'•12*57.2958 

00340  PRINT  "(INITIAL)  RADIUS  OP  CURVATURE  F:  I  =MR1 

00345  PRINT 

00350  IF  E<1  THEN  505 

00355  IF  EM  THEN  375 

00360  IF  2=1  THEN  375 

00365  GOSUB  765 

00370  GOSUB  510 

00375  IF  E  >=  2  THEN  560 

00380  PRINT  "ENTER  GANMA-H"; 

00385  INPUT  A 

00390  PRINT  "ENTER  HF-ii" ; 

00395  INPUT  D 
00400  PRINT 

00405  PRINT  "INPUT  FOR  HOBEL  PLUME" 

00410  PRINT  "GAHNA-M=" A , "HF-M="B 
00415  F=( (A-1 )/(A+1 )  >  0.5 
00420  Q=T2 
00425  R2=R1 

00430  PRINT  "FIRST  TRIAL  VALUE  ML  -1-1 " 

00435  PRINT  "ENTER  ML-M-1"; 

0 440  INPUT  3 
00445  I  =  FNC < B ) 

00450  .  =  P  *'  C  [it 
00453  (  =  ! 

•04  i  0  ETS'J?  2' 5 
00465  C3=: 

X-=J 

•')04?5  GOSUB  2  5 
00480  04=L 
O0->85  C -3-04+03 

00490  PRINT  "HL-M-1 ="8 , "THE TA-L-M-1  =  "C  *57.2953 

00495  E=E+1 

00500  GOTO  140 

00505  STOP 

00510  PRINT 


00515  PRINT  "PLUHE  SHAPE  CALCULATED" 

00520  PRINT  "X'V'R", "THETA" 

00525  FOR  S=0  TO  A  .2  STEP  0.2 
00530  T=FNA(SIN(T2)-(S/R1 ) ) 

00535  S=R1*(SIN(T2)-SINCT) ) 

00540  R=1+R1*<COS(T)-COS(T2) ) 

00545  PRINT  S,R,T*57.2959 

00550  NEXT  S 

00555  RETURN 

00560  IF  E>2  THEN  635 

00565  B5=B 

00570  R5=R1 

00575  PRINT  "ENTER  HL-M-2"; 

00580  INPUT  B 
00585  I=FNC(B) 

00590  X=I 
00595  60SUB  215 
00600  03=L 
00605  C - Q — 04+03 

006 1 0  PRINT  "SECOND  TRIAL  -'ALUE  ML  -N-2" 

00615  PRINT  "ML-N-2="B, ""HETA-L -N-2' "C*57. 

00620  E-E+1 

00625  B6=B 

00630  GOTO  140 

00635  B6  =  B 

00640  R 6  =  R 1 

00645  R6-R2 ) * ( B6-B5) / ( R5-R6 ) 

00650  I=FNC(B> 

00655  X=I 
00660  GOSUB  215 
00665  03=1 
00670  C=Q-04+O3 
00675  PRINT 
00680  PRINT 

00685  PRINT  "ML-H="rV‘THETA-L-M="C*57.2958 
00690  IF  APS < B-B6 )  0.001  THEN  730 
00695  B5=B6 
00700  R5=R6 
00705  GOTO  140 

00710  PRINT  "INPUT  THETA-F,R1 
00715  INPUT  T2,RI 
717  T2=T2/57.296 
00720  Z=E=1 

00725  PRINT  "PLUHE  GEOMETRY  SPECIFIED" 

0073 0  PRINT  "THETA-F="T2*57.29 58, "R1="R1 

00735  PRINT 

00740  GOTO  350 

00745  2=0 

00750  GOSUB  765 

O0,55  GOSUB  510 

00760  GOTn  -on 
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00765  P=n-J*F*2)'l(A/(A-t)) 

00770  PRINT 

00775  PRINT  "PRESSURE  RATIO  P=“P 
00780  PRINT 
00785  RETURN 

00790  PRINT  "THIS  IS  THE  FINAL  RESULT" 
00795  PRINT 
00800  END 
READY. 


Program  GRTRAN1 


GRTRAN1  is  a  laminar-transitional -turbulent  boundary  layer  program 
based  on  the  methods  of  Gruschwitz  [25]  and  Culick  and  Hill  [24], 


input 

quantity 

expl anation 

Y 

GO 

specific  heat  ratio 

Pr 

P 

Prandtl  number 

Mref 

MO 

reference  Mach  number 

p0 

PU 

stagnation  pressure  (psia) 

T0 

x/r*  vs.  M 

TO 

stagnation  temperature  (R) 

Mach  number  distribution 

in  the  nozzle 

output 

X/r*.  M,  9, 

HC,  6* ,  RI 

along  the  nozzle  wal  1 

Re,6  (test 

section),  5, 

6*  and  9  at  the  nozzle  exit. 

where 

HC  -  compressible  shape  factor  5*/9 

RI  -  equivalent  incompressible  boundary  layer 

Reynolds  number  (based  on  displacement  thickness) 


5  fs-INT  "PROGRAM  6»TRANt.6AS.  MHK  11/20' 1983" 

.0  PRINT  "TRANSITIONAL  BOUNDARY  LAYER,  GPUSCHWITZ-CtJI.IOK/H ILL-HHK. 
i?  PRINT  "MERGED  PROGRAMS  GROT  AND  B(.T2,  HHK  10/28/33" 

TO  00=3.141592 

25  DIM  X ( 50  > . M ( 50 ) . U (50 ) . D ( 50  ■ . S f  50  1 . K  50 1  .  H<  50) ..  B  <<  U  A  ■  '  0 
TO  D 1 1"  R '  50 ) ,  0  ( 50  i .  L  ( 50 ) ,  F '  50  * .  G  <  50  > .  n  f  50  * 

55  ofrtm  "L  AM  I  NAP  COMPRESSIBLE  BOU^D-Pf  LA'S8’  1  GRL'SLNUI  fZ)* 

TO  DEF  f f‘ 3  <  G  1  =  2 1 f n «  g  i  *  ( 1  - o ; *  5 r.  2 / 2 4 0  • 

15  pEF  FN5 ( 5  • =PNI  ( G )  -P* ( 1 2+5  "  2 ' 2 1 60 
:5 0  DEF  FNC  ( G  *  =  1  -  <  GO-1  1 F M p ;  r, .  /  *  j  l:  1 " ‘< F>  •  r>  1  1 

",5  DEF  P  *J  D  t  f*  >  r  1  fi  0  **  1  >  T  •  <  -  P  N  ft  '  <  -  '  /  1  2  ■ Nh  *  tji  <  <  l 

s 0  DEF  FNE  <  G  >  =  <  G 0 -  1  /  +  < 1  -  *  - *-i  7  r  <3  •  -  •  >  2  <•  ft  ■ " /  .2 I  -» 0  1  /  <  2  •* F N w «  ft  •  f  * 

55  DEE  F  N  T  I  ll  j  s  t  +  i  |,  0  -  1  i  !■:  ft  4:  rt  .  7 

>0  DEF  FNX <  Z  1  =Z 3’ /-3i  5-Z  /  045-Z  ‘Z  2ii’2  ■  2 

.'5  DEF  F  W  M  <  F  i  z  A  A  0  i  F  '  2  *■  N  0  ^  **  r.  -  .  i  -  i  u  - 1  it  H  iV ; .  >1 1" '.  3  I  1  -  ft'  2  !  ’■ 

10  DEF  P  n  N  <  M  ) =  <  F  N  T  <  M  0 )  /  F  -j  T  <  *1  <  <  ■  0  - 1  '  i  G  0  - i  '*  ' 

■(  S  n  P  P  P  M  r,  <  ”  )  r. '  i  •  * .  1  c  - '  ‘ .  v  j  T  -  :  ■ 


•i / ;;  i  .  9 


*0  AFP  FNh  ■ 


=  1?  3 1  5-ft  •’15-p-‘-j  /pm  • ; 

;  ,  '  o  :  •  3  ft  -34ft,  .713  !:• 


5  •  !j  >  )  t*  1  •  M  » rt  »  « ! :  '  •  .  )  : 

’  I  ft  A  1*  i  ■> ■>, 
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i  '(i  -  i3  T  *1 7  "EN'FP  PFANI'TI. 

i  -j  5  r  <j  p  tj  t  p 

!  PRINT  ''FfjTpp  MO  >=PFC.  PI  AIL  H  *  .  -  * 

I  35  :T  70 
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5 A  ’  F  ft.  f  '•  "  V  ''  1  P  V  ’25 

1 55  point  "E'jT-ftft  PE2D-'  •.,.3*1.  'mh,0  •  'IP  }  ll-  •LO.vjri.'D  '  .1  ;  ta  ;,  c  • 
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;  2  1  ft  ij  t  21  235 

>25  print  "ENTER  RO  ■’  =>i.  '■  V-' .  *:L  « 1  j . I ■  Ei:  • '  : 

230  !Nft"j-  v  » 

'  3*i  0  = .  ’  p 

’40  POINT  "I'-Annas  ’  A  Mll'iT-  ;  =  •  >  '  ■'  i«S;  ifti  3-'i  I 
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250  PRINT  “INPUT  MACHN  (YES), OR  U/UO(NO)=": 

255  INPUT  B$ 

240  PRINT  "ENTER  0  FOR  STAGN.  POINT.  1  FOR  SHARP  L.E."; 

245  INPUT  Z 

270  PRINT  "INPUT  N  OF  POINT  PAIRS="J 
275  INPUT  NO 

280  PRINT  “NUMBER  OF  POINT  PAIRS  X.M  (OR  U/U)="KO 

281  PRINT  "INPUT  FROM  FILE  GBATA  < Y/N )=“ ; 

282  INPUT  C$ 

283  IF  C$<>"N"  THEN  1441 

285  FOR  1=1  TO  NO 

286  IF  B$<>"YES"  THEN  297 

287  PRINT  “I="I ; “ENTER  X.M"; 

290  INPUT  X(I>.M(I) 

292  A< 1 , 1 ) =X( I ) 

295  A < 1 . 2 ) =M< I ) 

296  GOTO  300 

297  PRINT  "POINT  "  I"  ENTER  X  ( I  >  ,IJ/U(  I  >  =  “  J 

298  input  xm.um 

299  M ( I  )=FNM(U( I ) ) 

300  A( I , 1 )=X( I ) 

301  A( 1 , 2) =H( I ) 

.304  PRINT  ,,I  =  "I;,,X=,,X<I>:"N«"A(I.2> 

305  NEXT  I 

306  PRINT  "SHALL  BATA  A < 1 . [ ) , A ( 2 ,  T '  BE  STORED  IN  FILE  5DATA  (Y/N)?" 

307  INPUT  Bl 

308  IF  Bt="Y"  THEN  1470 
310  IF  fU  1)00  THEN  580 
315  IF  ZOO  THEN  580 

320  X< 1  )=  N< 1 ) =  Dll):  FT  =0 
330  PRINT  "NEAR  STAGNATION  POINT" 

335  FOR  1=1  TO  3 

340  PRINT  "I="I,"X=,,X(I);"H="M(I) 

345  NEXT  I 
350  F2=FNL(N(2) ) 

355  K ( 1 )=7. 700001E-02 
360  K=K ( 1 ) 

345  T0=7. 700001 E-02* ( X  <  2 )  -  X< 1 ) ) *FNN<N( 1 ) )/(R0+F2) 

370  D( 1  )=SQR( TO  > 

375  GOSUB  620 
380  B1=FNH(G) 

385  B2=FNI(G) 

390  PRINT  "THETA-0=“SRR( TO) ;"K="K; "ST AGNATION  POINT" 

395  X 3 = X ( 3 ) 

400  X 2 = X ( 2 ) 

405  F3=FNL(H(3) ) 

410  B( 2 )=B(  1  ) * (  1  - . 424*( ( F3/F2- 1 ) /( X( 3)/X (2) -1 > -1 ) ) 

415  PRINT  "THETA-2="D(2):"AT  X/L="X ( 2) 

420  H=FNM(F2) 

425  K2=FNJ<G) * ( BC 2) ) '2*R0*( F3 -F2)/(FNN( li( 2 ) )* (X( 3 )-X ( 2  » ) ) 

430  K(2)=K2 


435  1=3 
440  GOSUB  715 

445  1=4 

446  PRINT  "  X  N  THETA/L  HC  DELTA1/L  RI" 

450  GOTO  875 

455  D1=FNH(G) 

460  F2=FNL(M2) 

465  H=(M*M*(G0-1 )*FNB(G)/2+FNJ(G>*< . 3-G/60 ) ) /FMH(G) 

470  V< I )  =  1  +  < GO-1 ) *A< 1 ,2) "2/2 
475  H(I)=H 

480  F(I)=<H< I ) + 1 )/V<I)-1 
485  G( I  )=D( I ) /V( I ) "3 
490  PRINT 

492  M2=FNM(F2) 

493  D2=SQR( T2) 

494  D?=D2/D 1 
496  H2=H 

525  L(  I  )=R1*G(I  )"'2*<A(  1,2)  -A<  1-1 ,2)  )*V(  I )  ‘  .5/(A(  1,1  )-A(  I -1 , 1 ) ) 

530  IF  L<I)>=,02  THEN  545 

535  N( I )=1 O' ( 1 . 1 94*SIN( (L( I )-.02)*4.84*Q0) )*1 550/F ( I ) 

540  GOTO  550 

545  N ( I)  =  1 0 * ( 1 . 036 *S INI  (L <  I > - . 02 *6 . 06*QO > >*1 550/F ( I ) 

550  R( I )=R1 *G( I )*A( I ,2) 

551  T1 =R( I ) *V( I ) " (-3.5) 

555  IF  R( I )*V( I ) " (-3. 5 ) >N( I )  THEN  985 

560  PRINT  "T,U/T ,0="FN J(G)*( 1 ♦< GO-1 )*M0“2/2)/( 1 +(60-1 >+M‘2/2) 

561  print  x2;m2;B2;h2;d?;ti 
565  IF  G< - 1 2  THEN  595 

570  PRINT 
575  RETURN 

580  PRINT  "SHARP  LEADING  EDGE" 

585  PRINT  "Ml  = "HO; "AT  X1="0 
590  GOTO  860 

595  PRINT  "LAMINAR  SEPARATION" 

596  PRINT  "TURBULENT  REATTACHHENT  CONSIDERED" 

597  R( I ) =RO*G( I )*A( I ,2) 

598  T=I 

599  GOTO  1010 

600  STOP 

605  DIM  Z C 200 ) 

610  PRINT  "ENTER  K="; 

615  INPUT  K 

620  IF  K<?.479999E-02  THEN  635 

630  K=9. 469999E-02 

635  IF  «>-. 157  THEN  645 

640  K=-.157 

645  Z0= . OOOOOI 

650  TO=FNX(ZO)-N 

655  IF  K<0  THEN  670 

660  Z1 =1 2 

665  GOTO  675 


470  Z 1 =— 1 2 

475  Y1 *FNX(Z1 )-K 

480  Z2=ZO-YO*<Z1 -Z0 )/( Y1-YO) 

485  IF  ABS < Z2-Z 1 X.01  THEN  700 

490  Z1  =Z2 

495  GOTO  675 

700  G=Z2 

705  RETURN 

710  STOP 

715  XI =X< 1-1  ) 

720  X2=X  ( I ) 

725  T1»(D(I-1))2 
730  Ml I- T  > 

735  N2=H  ( I ) 

740  U1=FNL(«1) 

745  U2=FNL < M2) 

750  U9»(U1+U2)/2 
755  M9= 'HI  +f12>/2 
760  K1  =K(  1-1 ) 

765  B1  =  1 
770  K*K1 
775  GOSUB  420 

780  T2=T1+2*1X2-X1 >*FNN(M9) MFNA(G) -K*(2-rt9*H9*FNC<G) )/B1 )/<R0*FNL<«9) ) 
785  K2=81*R0*(T1+T2)*(U2-U1  ' /  <  2 *( X2-X! )*FNN<N 9)) 

790  K=K2 

795  IF  K2<- . 1 567  THEN  595 
GOO  GOSUB  620 
GO 5  N=N9 
810  B2=FNJ(G> 

815  IF  ABS < K2-K 1  X.0001  THEN  835 

820  B1  =B2 

825  K1=(K2+K1 )/2 

1330  GOTO  770 

835  D(I)-S0R(T2) 

840  K ( I ) =K2 
845  PRINT 
850  RETURN 
855  STOP 
840  1-2 

845  EK1  )=  K ( 1 )  =0 
867  X  ( 1 )  =  A  (  1 , 1 ) 

870  M< 1 ) =«0 
875  PRINT 

885  X( I )"A < 1 ,1 ) 

886  rt ( I ) =A( 1,2) 

910  GOSUB  715 
915  D  ( I  —  1  >=D(I) 

920  GOSUB  455 
925  K( 1-1 )*K( I ) 

930  1=1+1 

935  IF  I>NO  THEN  '780 


940  GOTO  875 
980  STOP 

985  PRINT"PQINT  OF  INSTABILITY,  TRANSITION  STARTS“I 
990  T*I 

995  D8=1 .43-1 /< . 02086* ( L OG < R < T ) ) ) *2.485) 

1000  H8=F(T)-D8 

1905  L8=(H8-.045) “ (-5.556)-. 185 
19 JO  T(6)=A( T,2) 

1015  T (7)=A(T, 1 ) 

1020  T(8)=B(T)*L 
1924  T ( 1 )=I 

1025  T( 2  >  =  N0 

1026  T(3)=Rt 

1027  T(4)-H8 
1928  T(5)=L8 

1029  T (9>=L 

1030  PRINT  "TURBULENT  COMPRESSIBLE  BOUNDARY  LAYER,  AIR,  HHK ,  8/24/81" 
1035  PRINT 

1038  N=4 
1040  RO  =  T ( 3 ) 

1045  I'  =  T(  9 ) 

1050  M0  =  T<6) 

1  055  XO=T<  7) 

1  060  NO  =  T  < 2 ) 

1065  D0=T(8) 

1070  IF  DOOO  THEN  1080 
1075  00=  9.999999E-IO 

1080  PRINT  "EXPONENT  OF  SHEAR  LAY,  N* ( TENTATIVELY) ="N 
1  090  C0=(N+1  )*.1266*N'(-2.6815)/(R0'(1/N)) 

1  095  S0=(  (D0/D)*(H0/( 1  +  .2*H0*N0) ) "3) " ( 1 +1 /N) 

MOO  PRINT  "TRANSITION  CONDITIONS  XO(FT)=MX0;"M0-’,M0;“''0="D0 
1105  PRINT  "SHEAR  LAU  EXPONENT  N="N:"REY-0,D="R0 
1 1 07  PRINT  "  X  N  THETA/L  HC  DISF7L  RI" 

1110  LO  =  T ( 5 ) 

1115  PRINT 
1120  X1=X0 
1 1 25  (11=110 
1130  D 1  =  D 0 
1135  J1=0 
1140  11=0 

1145  Q0=( ( (D1*( 1 ♦ . 2 +HO*MO ) * (-3 ) /D ) "  ( 1 +1 /N ) ) /CO i *4 
1150  Q1 =Q0 
1155  Z1 =0 

1160  FOR  0=T( 1 ) + 1  TO  NO 
1165  X2=A(0, 1 ) 

1170  H2=A<0f 2) 

1175  IF  X2=0  AND  N2=0  THEN  1299 
1180  H9=(H1*N2)/2 

1185  12=  ( X2-X1  )*(M?"<3+2/N)*n  +  .2*H9*M9)M(1/N)-4>) 

1190  11=11+12 

1 195  B2=B*< 1 ♦ . 2*N2*N2)  '3*< (C0*1 1 +S0 )*H2  " (-3-3/N) ) * (N/(N+1 » ) 


1200  D9=(D1+D2)/2 

1205  Q2=m/M0*(-3*(1+1/N))+Q0  (1/4))  4 
1210  R2=R0*D9*M9/((<1+.2*M9*M9) *3>*D> 

1215  B=.0304*L06(R2>-. 23 
1220  J2=(B-L0G<M9/M0))*(Q2 -01 ) 

1225  L2=Q0*L0/Q2+L0G(H2/»O)+(1/«2)*(J2*,Ji ) 

1230  IF  L2<-. 1 5  THEN  1400 
1235  IF  12 <-.185  THEN  1410 
1240  K2=.045+(L2+.185)A(-.18) 

1245  H2=(K2+1 ) * ( 1 +. 2*M2*M2 > -1 
1  250  T1=R9*D2*M2/<  ( <1*.2*H2*M2)A3)*D) 

1255  0°=H2*D2 

1260  print  x2:H2;D2;h2;d?:ti 

1265  XI =X2 

I2?0  Ml =  M2 

1275  Jl = J1 +J2 

1280  CM  =Q2 

1285  NEXT  0 

1  290  T1  -RO*D2*M2/(  •'  ( W. 2-*M2*r'2 )  " 3 >  k D  > 

1295  PRINT  "  R ,  T  =  TEST  REYNOLDS  t»="T1 

1300  PRINT  "SELECTION  OF  SHEAR  LAW  EXPONENT  'HERE  USED  N="N" > " 

1  305  PRINT  "POP  1 00<R , T < 5 E 3 ,  N  =  4" 

(310  PRINT  "FOR  500  R.TC5E4.  N  =  5" 

1315  PRINT  "POR  3E3<R,T<6E5,  N  =  6" 

1320  PRINT  "RECALCULATE  UITH  PROPER  4  IE  NEEDED" 

1325  N1=2*N-1 

1330  C2=M2*M2/(5+M2*M2> 

1335  Y=1 

1340  Y0= 1 / ( N1 + 1 ) - 1 / (N1 +2 ) 

1345  Y  =  Y  +  1 

1 350  Y2=C2  ■ C Y-1 )/(N1 -1 +2*Y)-C2* ( Y-1 )/(N1+2*Y) 

1355  Y0= Y0+Y2 

1360  IF  Y2  .0091  THEN  1370 

1365  GOTO  1345 

1  370  03= ( 1 -C2  >*N1 *Y0 

1375  PRINT  " D E I. T A/ T HE T A = "  1  / 0 3 

1380  PRINT  "ROUNDARY  LAYER  THICKNESS  (FT»-"02/D3 

'385  PRINT  "BOUNDARY  LAYER  MOMENTUM  THICKNESS  (ET)="D2 

1390  PRINT  "BOUNDARY  LAYER  DISPLACEMENT  THICKNESS  (FT)s"52«H2 

1395  GOTO  1415 

1400  PRINT  "SEPARATION  IS  POSSIBLE" 

1405  GOTO  1235 

1410  PRINT  "SEPARATION  MUST  HAVE  OCCURRED" 

1412  G0T01325 

1415  PRINT  "IS  NEU  SHEAR  LAW  EXPONENT  DESIRED  FOR  RECALCULATION?! Y/N)" 
1420  INPUT  B$ 

1425  IF  8$<>"Y"  THEN  1540 

1430  PRINT  "ENTER  NEU  SHEAR  LAW  EXPONENT" ; 

1435  INPUT  N 

1440  GOTO  1040 

1441  FRENI ="GDATA" 


1442  FOR  1=1  TO  NO 

1443  INPUT  II,  A<I,1),A<I.2) 

1444  X ( I ) -A ( I  •  1  > 

1445  H< I )=A( I ,2) 

1446  NEXT  I 

1447  CLOSE  N1 

1448  FOR  1=1  TO  WO 

1  449  PRINT  “T  =  ,,i:"X="*(I,1);,,H=-A(l 

1450  X(I)=A(I,1> 

1451  NEXT  I 

1452  GOTO  310 

1  470  FILEN1  ="ODATA“ 

1480  FOR  1=1  TO  NO 

!  490  PRINT  HI,  IJSINO'HH  .4H4H1' ;  A  ( I  „  1 

1510  NEXT  I 

1520  CLOSE  41 

1530  GOTO  310 

1 5^0  END 


Program  BARTZ 


BARTZ  is  a 

fully  turbulent 

boundary  layer  program  applying  the 

Bartz  method  [18] 

• 

input 

quantity 

symbol 

explanation 

Y 

K 

specific  heat  ratio 

R 

R8 

gas  constant  (ft  1 b^/ 1 bm  R) 

p0 

PO 

stagnation  pressure  (psia) 

T0 

TO 

stagnation  temperature  (R) 

u 

MO 

stagnation  dynamic  viscosity 

(lbm/sec  in.) 

r* 

RO 

throat  radius  (-) 

R* 

KO 

throat  radius  of  curvature  (-) 

Cp 

Cl 

specific  heat  at  constant  pressure 

(Btu/1 bm  R) 

k 

K1 

wall  conductivity  (Rtu/in.  sec  R) 

L 

wall  thickness  (in.) 

h 

HI 

film  coefficient-coolant  side 

T9 

coolant  temperature  (R) 

x/r*  vs.  r/r* 

nozzle  wall  contour 

output 

x/r*,  r/r*,  5**/r*,  6/r*,  6*/r*,  DELT/6,  TW/TO,  HG  along  the  nozzle 
wal  1 


at  the  nozzle  exit 


~  ? A?fTZ -r’JLLY  TlJWJlzNT  BOUNDARY  ‘_4'?CC" 

10  ::«  X ( 25 , 3 ) , Af  <  3 ' , B  S ( 3 ) , A ( 2 5 . 7 1 
20:  88.888  88.888  ##.##8  88.888  88.88#  '‘*.8*8  88.888  88.8^8 
30  DEF  FNX(M>=<2/(K+1)+(K-1)*M*M/<K+1)>A<(K+1)/<2*(K-1)))/M 

40  DEF  FMY (H)-(K+1 )*M/(H*N*(K-1 )+2)-1/M 
50  DEF  FNA(H)*T*(1+(K-1 )*H*M/2) 

AO  DEF  FNB(M)=<1/T-1)*VA<-1/7) 

70  DEF  FNC(H)=(K-1  )*H*M/ ( (2+<K-1 ) ) *T ) 

80  DEF  FND(H)=FNC(M)*MA<2/7) 

90  DEF  FNStH)=1/(<0.5*T*<1+<K-1)*H*M/2>+0.5)-(O.A)*(1*<K-1 ) *M*M/2) A (0. 1 5) ) 
100  PRINT  "TURBULENT  BOUNDART  LAYER  IN  COOLED  C-D  NOZZLE" 

110  PRINT  “BARTZ  METHOD,  MODIFIED  BY  HHK  7/1/76" 

120  PRINT  "ALL  BL.  THICKNESSES  DIMENSIONLESS  IN  TERMS  OF  R-THRQAT" 

130  REM:  NOTE  THAT  INPUT  DATA  ARE  IN  UNITS  CONFORMING  WITH  BARTZ 
150  PRINT  "ENTER  GAMMA, GAS  CONST,PO,TO,MO,RO, KO"; 

160  INPUT  K,R8,P0,T0,M0,R0,K0 
170  PRINT  "ENTER  C-P,K1 ,L,H1 , T9" ; 

180  INPUT  Cl ,K1 ,L,H1 ,T9 

190  PRINT  "ENTER  NOZZLE  GEOMETRY  X/RO,R/RO" 

200  PRINT  "NOTE  THAT  THROAT  <R/R0=1)  MUST  BE  ONE  OF  LOCATIONS" 

210  :88.  8888888.888  888888.888 
220  PRINT  "ENTER  8  OF  POINTS  (<25)"| 

230  INPUT  N 

240  PRINT  "NUMBER  OF  POINTS="N 
250  PRINT  "P0INT8  X/RO  R/RO“ 

240  PRINT  "GEOMETRY  FROM  FILE  3?“; 

270  INPUT  A* 

280  IF  AI="YES"  THEN  2250 
390  FOR  1=1  TO  N 

310  PRINT  "ENTER  X( "I") , t > ,X( "I" ,2' 

320  X( I , 3 ) =1 

340  INPUT  X(I,1),X(!,2) 

330  PRINT  USING  20, I ,X < I , 1 ) ,X< I ,2) 

340  NEXT  I 

370  PRINT  "SHALL  DATA  BE  STORED  IN  FILE  3?"; 

380  INPUT  B$ 

390  IF  Bt="YES"  THEN  2300 

400  PRINT  "GAMMA="K;"GAS  CONST. ="R8 

410  PRINT  "STAGN.  PRESS. ="P0;"STA5N.  TEMP="TO; "STAGN .  VISC.="MO 
420  PRINT  "THROAT  RADIUS="RO; "THROAT  CURV.="KO 
430  PRINT  "C-P="C1;"K="K1;"L="L;MH-E="H1;"T-E="T9 
440  U0*(2/( K+1  ))A((K+1)/(2*(K-1)) )*PO*SQR (K*32. 1 74/ (R8*T0 ) ) 

430  R9=R0*U0/N0 
440  P=4*K/(9*K-5) 

470  N1»C1*M0/<R0*H1) 

480  N2=L*C1 *M0/(R0*K1  ) 

490  N3=0. 0228*R9* ( 3/4 )/P" (0.44) 

500  PRINT  "REF.  REYNOLDS  »="R9; “PRAMDTL  #="P 


510  VI =PA ( -0. 34 ) 


[•r.i 


07 


520  GOTO  1130 

530  A0=R*2 

540  IF  C=0  THEN  660 

530  H8=1 .4 

5 60  GOSUB  580 

570  GOTO  630 

580  A1=FNX(H8) 

590  IF  ABS < A0-A1 XO.OOOI  THEN  630 
600  M9=H8+(A0-A1 )/( A0*FNY(N8) ) 

610  H8=N9 

620  GOTO  580 

630  M=H9 

640  RETURN 

650  STOP 

660  M8=0.01 

670  GOSUB  580 

680  GOTO  640 

690  11=0 

700  Z1=0 

710  00=0.05 

720  Z2=Z1 +D0 

730  Z9=Z1 +00/2 

740  IF  Z 9> 1  THEN  780 

750  11  =  11  +Z9A7*  ( 1-Z9)  *=D0/  ' 1  +B0*Z9-C0*Z9*Z9 ) 

760  Z1=Z2 

770  GOTO  720 

780  RETURN 

790  12=0 

800  Z1=0 

810  00=0.05 

820  Z2=Z1 +D0 

830  Z9=Z1+D0/2 

840  IF  Z 9 > 1  THEN  880 

850  I2=I2+Z9“7*B0/(1+B0*Z9-C0*Z9*Z9) 

860  Z 1 =  22 
770  SOTO  820 
880  RETURN 

890  IF  ?*!<•  i  THEN  930 
00  F  =  ?  -1 ,25*(  1+S1 ) •* ( X2— X 1  )/SQR(  < K + 1 ) *K0/2) 

910  N0=I 
920  GOTO  940 

930  F=1-5*(M1*H1-2*S1-3)*L0G(R2/R1 ) /( 4*( 1 -Ml *H1 ! ) 
940  F0=0.0285*FNS<H1 ) * < W 1 *R 1 *R 1 /R9 ) ' ( t /4) * <X2-X1 ) 
950  F2=(Fri.25*F+F0)'(0.8) 

940  RETURN 
970  B0=FNB(N) 

980  C0=FNC(N) 

990  RETURN 
1000  B0=1/T-1 
1010  CO=FND(H) 

1020  RETURN 


1030  G 1 =F 1 *K8*( 1 -T7)/ (I8*IM ) 

1040  G2=F2*K9*( 1 -T )/( I9+R2) 

1050  S0=1 +?*(L0G<G1 ) -LOG ( G2 ) )/8 

1060  G3=!»,02565*FNS(H1 )*I8* ( W 1 /R? ) A0.25*SGR ( R1 )*(X2-X1 )/(K3*F1  '  1 .25*P* ( 0.46) ) 
1070  G4=G8*(9/7)*G0+G3 
1080  G9=<ABS(G4) ) ' ( 7/9) 

1090  H4«H3*F4S  (H2)  *R2*  ( -3/2)  *«2 '  ( 1  •' 4 )  *G9 "  ‘  - 1  /? )  / F 2 ''  ( 1.'4 ) 

1100  T 2  = I  - i 1 -T9/T0)/( 1 +N4* (N1  +N2) ) 

1110  PRINT 
1120  RETURN 
1130  Q=1 
1140  C=0 
1142  Q=Q 


1150  FOR  1*2  TO  N 
1160  X 2= X (1,1) 

1170  R2=X<  1,2) 

1180  R=R2 

1190  IF  Q  2  THEN  1920 
1200  T=T7 
1210  T 1 =T7 
1220  GOSUB  1240 
.'230  GOTO  1540 
1240  R=R2 

1250  IF  R=1  THEN  1290 
1260  GOSUB  530 
1270  M2=M 
1280  GOTO  1300 
1290  M2=1 
1300  V=G8 
1310  GOSUB  970 
1320  GOSUB  690 
1330  19=11 
1340  GOSUB  790 
1350  J9=I2 
1160  COSUB  1000 
1370  GOSUB  690 
1380  K9=I1 
1390  GOSUB  790 
MOO  L9=I2 
1410  U2=7*I9/FNA(rt2) 

1420  S2=(FNA(N2) /7-J9)/I9 

1430  GOSUB  890 

1440  IF  Q=2  THEN  1460 

1450  GOSUB  1030 

1  460  IF  0  13  THEN  2<30 

1  470  RETURN 

1480  GOSUB  ■ 730 

1 4°0  IF  R2M.00001  THEN  1520 

1500  C=1 

1510  Q=Q 

1520  NEXT  I 


1530  GOTO  2220 

1540  IF  Q<3  THEN  1710 

1550  IF  A8S(T2-T1 ) > 0 . 0 1  THEN  1880 

1560  A(I,1 )=F2 

1570  All, 2)=F2*FNA(N2 )/( 7*19) 

1580  A(I,3)=S2*F2 
1590  A( I, 4)=G9 
1600  A(I,5)=T2 
1610  A(I,6)=N4*CI*MO/RO 
1420  A(  I,7)=rt2 

1430  PRINT  "  X/RO  R/RO  DEL*  *7R0  DEL/RO  DEL */R0  DELT/DEL  TU/TO  H-G" 
1640  PRINT  X2;R2;F2;F2*FNA(H2)/(7*I9);S2*F2;G9;T2;N4*C1*N0/R0 
1660  PRINT  "COMPLETES  ITERATION  FOR  POINT"G 
: 680  PRINT  "UHERE  R/R0="R2; "AND  MACH»="M2 
1490  T8=T2 
1700  GOTO  1 "20 
■<0  T8-T2=T7 
•>  ’20  GOTC  14CJ 
|  770  *  1  * <2 


!  740 

PI  c  P  2 

1  750 

o=n*i 

1760 

FI  =F2 

1770 

G9=G9 

1780 

Ml  =M2 

1  700 

18  =  19 

1800 

J8-J9 

1810 

K3-K9 

1820 

L  8=L  9 

1830 

IF  G<3  THEN  2140 

.'840 

T7=TB 

1850 

SI  =S2 

1860 

U1  =U2 

1870 

RETURN 

i  880 

T1=2*T1 *T2/ ( T 1 +T2 ) 

1890 

T  =  T 1 

1900 

GOSUP  1240 

1910 

GOTO  1540 

1920 

FI  =0 

1930 

V*V1 

1935 

X2=X( 1,1 ) 

1936 

R2=X ( 1,2) 

1940 

GOSUB  530 

1950 

Ml  =H 

1960 

G8=V1 

1970 

T=T7=0.25 

1980 

GOSUB  970 

1990 

GOSUB  690 

2000 

18=11 

2010 

GOSUB  790 

2020 

J8=I2 

7030 

GOSUB  1000 

GOSUB  690 
K8=I1 
GOSUB  790 
L8=I2 

W!=7*I8/FHA( HI ) 

S1=(FNA(M1 )/7-J8)/I8 

QSQ*1 

X 1  =X2 

R 1  =  R  2 

GOTO  1150 

G?=r 

GOTO  10*0 

T7  =  T7 

GOTO  1900 
G?  =  G3 

N4=0 

PRINT  "POINTS  1X2  GENERATE  INPUT,  EXEPT  FOR  PEL**/RO  AT  POINT 
GOTO  1 470 

R8=R?*F2*FNA(N2)*'H(K-1 >t«2*H2/2)  0.6/<7»I9*R2*R2> 

PRINT  "AT  LAST  STATION,  CEYNOLDSN  RE, FREE  STREArt, DELTA="R8 

GOTO  2340 

FILE  HI s"BDATA" 

FOR  1=1  TO  N 

INPUT  t»1fX'I,3),X(I,n,X(I,2) 

NEXT  I 
GOTO  400 
FILEN1  ="BDATA" 

FOR  1=1  TO  N 

PRINT  Nl,X(I,3),Xn,1i,X(I,2) 

NEXT  I 
CLOSE  #1 

PRINT  "NOZZLE  COORDINATES  STORED  IN  FILE  BDATA" 

GOTO  400 
END 


‘0*4. 


Y 

G 

specific  heat  ratio 

n 

N 

flow  profile  exponent 

ml 

Ml 

exi>t  lip  Mach  number 

MJS 

M2 

jet  surface  Mach  number 

Re,6 

R1 

Reynolds  number  based  on  velocity 

boundary  layer  thickness 

output 

Shear  stress  and  eddy  viscosity  functions  and  spread  rate  parameters 
along  the  plume  boundary. 


5  PRINT  "PROGRAM  FREES-DEVELOPHENT  OP  FREE  SHEAR  LAYER" 

10  DIN  AT 5. J 00) 

15  P0=3 . 1 4 1 592 

20  PRINT  “THIS  VERSION  WITH  PLOT ,HHK  9 '30/80" 

30  PRINT  "LINEAR  NIXING  PROFILE  IN  EXPANDED  1 'N  PROFILE  D.L 
40  PRINT  "THIS  VERSION  OF  8/28/80  CAN  HANDLE  THIN  AF'PR.  8.L 
30  DEF  FNZ'P)=f  N/S3PM-.:  i  -cr'2*F  2)*P  ( (G-1  ' /G) ) 

60  DEF  FNC(rt  >  =SQR(N  2.  ‘ 2/ ( G-’ ' -H  2  >  '■ 

20  DEP  FHNf  P)=F  N*FHD(F ) /i 1 -F  tr  *C ’ *C1 ) 

80  DEF  FND(F)  =  ( 1  - < 1 -Cl  ' 2+F"2 ' *P '  • (O-1 '  -  G  • )  0.3 
RO  DEC  FNA(P  '  =c  lJ  ' -F  2*Cl  '2  • 

1 00  DEP  PNVfP  '=D2M  - F 3 L 0 G  L  f  1  *C0'/(  I  -  1 0  1 '  ■'<  2*C3> ) 

120  PRINT  "GANNA=": 

•  30  I«PtjT  G 

no  C'S'JWT  "PROFILE  EXPONENT  M  =  ": 

150  INP'L'T 
I  ,0  pRT»)T 

1 70  INPUT  0! 

I  30  PRIMP 
■|«0  INPUT  62 

200  F'PIMr  " 4  OF  INTEGRATION  STEPS-" : 

210  INP'JT  D 

220  c«  I«T  “Enppr  S'EvN2LD'F  «•  *  -  A-DEL  r4  ' 

,:30  i  M p ij t  ri 
24rt  K=i  2  +  2. 7t  - 'f  2 
250  2"syO  =  P8  =  .'s1 

240  R?=  1 0 " '  4 . *-0 . 039  *  A2 >-0 . 053*42  .2 - 
270  :■=  i  ( 1  *(G-i  '  'NT  2 -’2 )  / •  i  + 1  G -i  •  *M2  ' 2  '2 "■  1 
280  P&INT  "GAnNA»hGT"L'PSTREA«  5.L.  PROFILE  EXP.--“N 
2°0  PRINT  "N1  =  "N1 ; "N2S"N2 :  "P'2  F”  s"?‘ 

300  PP'INT  "REYNOLDS  *, ’■  1  -A-DE'_TA1  **"Ri;  "S  H?NA  •  CONST,  ■  =“K 
310  PRINT  "4  CF  INTEGRATION  $?EPS="B;"REYNOLDS  »,  Tt'AN'S.  ="R 
120  C !  =  r  ,J  C  ?  0  ’  ! 

■  jn  r-  2  s  f  r t  m  2  » 

’•40  J i  =  S -3 -  ’  - '  !  -  22  2  1  ♦EXP1 :  LOG 1  (  ■  *C2)  / ( 1 -C2) '  i  2-2  •  '■ / C 27 
JfO  PPI«T  "SIMILARITY  VALUE  FOR  J 0  = " J 5 
U 0  Z>  sPnR> •*>  1  -C  2'  '  2/C* 

!?:  » r.'srt 

590  » i  =  Z  3  = ' 

5°0  Z5=»£M0 
100  p  1  s F 

HO  IF  ASS^F1-F8) ■ .0001  THEM  -10 

420  D=D / 1 0 

440  C2=F 1+1/0 

450  F9=<Fi+F2'/2 

440  IF  P9 >*  THEN  ; 020 

470  A  =  P N [h  F 9  i 

480  12  =  11  ♦FNMi’F9)t.'F2-F’  • 

190  i'2  =  Tl  *p  -  ( -■  /«  •  *rJ».7 1  *Ffi’  e?if:F2-Fl  > 


# 

.\w:v 

^  ^  ;-v ?r yg ■.  rn^T^r^v.r*; 

i  ti- 

J  L 

5 

95 

5  , 

*  *  * 

*  .  * 

•500 

Y3  =  r2 

510 

Dt  =  It*P*(-t ) *W*C t /< -  f  +f_0G(  n+FNB  F  I )  >  /< 1 -FNt'IF; »  )  •  '  <  2 ■>  Fm<  *  1)  1 ) 

]  R 

S20 

D2= 1 2*P “  (-1  )*N*C!/<-1H06«  *  1  +FND(F2>  )/(  1  -FND ^ F 2 )  >  >  •'>  2*p«!0<F2  >  > ) 

J  B 

530 

IF  Dt  =0  THEN  550 

5  40 

Z7=2 >*C  1  »FHIi'-1  '  >N/IH 

;•  >; 

550 

1=2 *P " (-1 > *C’ *FMD(F2) *H'D2 

560 

J2  =  J!+FNA(F!»)*(F2-F1  ) 

570 

if  [n=o  twfn  5^o 

1 

580 

Z  8 = J 1  +  2  7 

!  P* 

590 

IF  P I =  0  THEN  680 

t  .'I- 

<600 

Q0=SSR <  1  - <  i  -  FNJ •'  F 1  1 )  '  2 '  >£*»'  78  ■  • 

fc.'  ' 

•  :'0 

P0=30/C2 

520 

j;= j1  +fnv  co  ,  t  •:  -  2-PI 1 

530 

P3=FNP(F1 1 

|  * 

.540 

e  j  i  n  o  /  0  3 

1 

6  50 

1 13  =  n  i  $  ( - F 4  >■  L O 6  ( <  l  *O0  > <  1  -Q>)>  :■  /  <  7  *  n  v>  % 

L  V 

b  ’/ 

560 

l'z jpti 

*;  '.- 

•'D  V 

T-  111  0  T^EH  .;?•} 

* 

U 

•880 

:«=o 

w 

s. 

5 ’0 

C  0  =  S  G  F  <  *  -  ■:  *  ~>FNDfF2’  1  I  >  *  t  *F'(  10 '  1 

1  f- 

'O.i 

3  -  C  0  C  2 

W 

MO 

;r  y.jsi  THEN  7:0 

4  «  •> 
4.  J 

-  PINT 

V  •:; 

’30 

print  "phi-i-'I  •s"Ft:"?m-i  :i*i  >  :''t-2 

«  '  4 

4 

•40 

-R  HT  "TMI  -  2>  I  >  =  "F  JO  •  - 1  - C 2 :  "•-■NT-  ><  r*i  '•s"FND(r2 VC2 

» 

'50 

=•*=■: »J T  H'2/Yt*» /::"PELTA/PELTA1="&2;P1 

in 

'  ,4  ij 

PRINT  "  J0="  JO  OO'.T  41  r  '  T  ;-p  2 

PRINT  J0-3ATIQ  =  ''J9/ J3 

’  5  "i 

1 3  =  ~  N  [i !  F  2 ) 

:;  >• 

’<>0 

IF  Z0  =  0  then  '450 

•  > 

:oo 

Zi-Cl* <  t- JO  2*C 2"2)  0  2 

■:  ■* 

510 

F  3  =  C  0  •'  C  3 

ip 

■'320 

>330 

Z2  =  t 75+18 : /2 

2 4  =  I> 2  •  1  OG ■  1 ♦  CO  VH  -CO >  "•  / ?2*C3  > ) 

>,  y 

■340 

T*  =§«•■.  e.,;  ,  +  •  C2  2 -  ♦  >  * •  Z4-Z«  ■>  ' 

• 

,->50 

E0-Tl  *  ! 1  -CO  2 ) * C 2 <■ 1 4*(  1  -C2*C2 '■ ) 

:  <• 

>35' 

IF  Z 0  = 1  t»£n  370 

►  i- 

!v) 

PRINT  "cucae;  STRESS  F U N C  T I O N 3 " T 1  ; C T .0 1’  -.'IS C 0 '3 1  r v  jnp ;  ihn 0 

t 

■J-.) 

V“sY1  f  .  74-7V.  tZI  fV> 

p  — 

»80 

GPS '33  1  ‘"’O 

o 

5?  00 

IP  TO-I  T^EN  «10 

■* 

?oo 

PRINT  "X'DELTPi,  l A « .  =  ' v 2 

•>  •  o 

R2  =  '  C2  pR1  *  <2  >P  ‘  ■'  M  -  Cl  2 1  -C7"  2  •  *  7)  -'Cl 

/  **\ 

■<:o 

PRINT  "PETNCLD?  P  ="R2 

f  ^ 

530 

IF  P 2  F:°  then  1 330 

► 

7  40 

23  =  24 

▼^1 


7.V." 


M 


'  000 

F!  =F2 

1010 

60  TO  410 

1020 

Y3=T2 

1025 

1= J 

i  030 

PRINT 

1040 

X1=X2 

1050 

^ PINT  'TOGE  r.F  RME4R 

Li-V'ER  RE4C4 

1055 

WO*  J 

'060 

fPINT 

i  }'0 

S  0  =  -  *  *L0G'  •  1  *C2  >  ; 

2  ■  ■  ";2‘C2) 

'090 

91 sN*Ci ‘12 '9-Y? » C 2  2 

i  1  .  r  2  2  > 

1090 

!  i  c  ;• 

T 1  =  C 1  M  ‘  J2  ••‘e- 1 2  »C2/  •.  1 

;2=  ;T 

-f*>*  ;) 

■mo 

rn=  <  S’  *92)  so 

•1  i  *{\ 

PRINT  -  Y  j  ' tl  T  il  =  “  r  ■? : 

"[i;-Lri  j'li  ‘ 

1 1 3  o 

t2  =  y  j‘C2  '  ‘ 1  -  C  2 "  2  '■ 

1  i.) 

•_  =  2»C2  m  T  !  -T? )  / n 2 

M  jO 

.0*SPR< *  -  ■ -c:  2  * k ■ : 

■  5  .  1  j 

1  1 30 

- PIMT 

J  1  ^  •*) 

»>r  i  i r j i  =  ■  •  i o 

> 

'  '90 

I  e  p  ;  ■■  p  o  t  c  fj  t  ,f  .  0 

i  ’Or) 

2;.  =  c:»  .!-"!•  1 0  ,."■ 

Tr> 

l_2_X0 

r  7  -  _  r  *  7 

;  If;1! 

2  2=  25* '*>  2 

j 

Z  a  =  ?  2  4 ....  r  :  ♦  ?  5  o .  ,  t  ..  r  •  • 

f  “• 

'210 

SMH!ZK:  >  ■  21  '22 

:  240 

R2='C2‘F  t  «*2*;,*«  1  1  -C1 

2  •«  ,  *-C?A^ 

•  130 

p-, . ..  v  . j_  ” 4 1 , 4 ** ,  = 

"  <2  : "  RE  *  "O'.: 

'  2  s  0 

I  2?a 

i  i  =  '<2 

•!  290 

Y  3  =  Y  3  »  2 

•’  2 "  0 

:p  9 1':  2E  ‘■06  THEN  io:o 

I  ; .)  m 

505 JB  t 7 90 

I3'0 

GOTO  '’00 

Uj} 

:  :■  ’  ,s  r 

i  330 

*  ~  Y  •)  •;  7 :  i |7  ^  !  d i)f 'i 

‘  1  ■*  0 

Fg=c: 

■t 

0 . 1 0  ...  1 . 

i  360 

i  7  7 1) 

Y  0=0  * 

'  390 

(\  —  ♦ 

'  ;*0 

GOTO  4 '!  (J 

1 400 

PRINT  "TP415  FTIOV  'J4S 

9 C CURED" 

M05 

Is  . 

1410 

10  =  0 

■i  420 

0  8  =  X  2 

1  430 

07  =  J 

;  440 

y  1  :  V  2 

1450 

'JOSHS  1  05:  . 

1  4*i) 

PRINT  "1>- lELTO' .  "MRS 

.11  V  ■> 

'^0 

T 1 =9  OR . °  )  ‘ 1 ’  / 22  '  3-  '  * 

, . :  j--;-  •. .  •  0 ' 

I  480  E0=T1  1-C0  2'  *22:  ( A*<  >  -L2‘C2.» ) 

I  490  PRIMT  "SHEAR  STRESS  FJNCT  [0rJ-,,T1 ;  "EOD  Y-VI'SCO'S!  T1  FUWP  T  £  ON  =  "  EO 
1500  GOSUB  1?’0 
1510  GOTO  960 


1520 

GOSUB  1550 

53C 

f.iSUB  1~c0 

5  i  o 

GQtq  1440 

i  550 

-  3  =  r o  'C3 

1  56  0 

Z  3  =  Z  ^ 

1  5’0 

!  530 

X2=<1  +  '  Z  »-:3  ■ 

!  5P0 

RE 'URN 

1  600 
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!  A 1  f) 

.  2 =  ( S  *  ♦  5  - *  'SO 

1420 
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Program  ASUEXP1 

Program  ASUEXP  calculates  the  sudden  expansion  and  impingement  of  a 
plume  on  a  solid  cylindrical  wall. 


input 


juantv 


re/r* 


rimping/re 


symbol 


explanation 


specific  heat  ratio 
nozzle  exit  radius 
nozzle  exit  lip  angle  (deg) 
nozzle  exit  Mach  number 
measure  of  the  local  rate  of 
acceleration  near  the  nozzle 
lip  [12] 

radius  of  the  impingement  wall 
Reynolds  number  based  on 
velocity  boundary  layer 
thickness 

velocity  boundary  layer 
thickness  at  the  nozzle  lip 


output 


S  £ 


Peak  pressure  ratio,  location  of  pressure  peak,  inviscid  impingement 
angle,  base  pressure  ratio  and  location  of  the  viscid  pressure  rise. 


This  program  was  written  in  HP-BASIC  and  run  on  a  HP  9830A. 


r.  r.  O 


100 


M 


1  con  CC  2?  3 

10  PRINT  "PLUME  IMPACT  ON  CYLINDRICAL  WALL'  ZERO  BLOW-E'Y.  THEORY  AND  E:  F . l i 
710  PRINT  "PLUME  IMPINGEMENT  ON  A  CYLINDRICAL  HALL.  HHl.S  31.7?" 

715  PRINT  "ITERATIVE  SOLUTION  POP  FINDING  ZEPO  BLOH-BY  CASE.  HHI  ?  5  7? 

730  PRINT  "CALCULATION  OF  APPROXIMATE  PLUME  SHAPE" 

735  PPINT  “FOP  IDEAL  OR  NON-IDEAL  CONICAL  C-D  NOZZLES" 

730  SIM  D»C31«  ICS. 31 
735  DEG 

740  OEF  FNA<::;.=ATN' SORM -:>:>1E-?S  ■  • 

745  Def  fnbo:..*atn.  scr<i-7»:  •  •  x+ie-s?  ■  ■♦3*atnie  ::  o* 

750  DEF  FNP' M  .*1  (3*.  A-1>+M*M>  .  A+l  )  t '  •  A+i  •  •  >..  A-l  •  •  •  M 
755  Z=1 
7 SO  V=1 

7C5  PPINT  "NOZZLE  FLOW  SPECIFIED" 

770  DISP  "ENTER  GAMMA*"! 

775  Input  a 

7 SO  PPINT  GAMMA* "A 

735  DISP  IS  NOZZLE  IDEAL  CONICAL'""'  5 

7 SO  INPUT  DI 

7S5  DISP  "ENTEP  NOZZLE  EXIT- THPOAT  RADIUS  PATIO*": 

Boo  input  ps 

S 05  DISP  ENTER  NOZZLE  EXIT  ANGLE*"  ’• 
sio  Input  c 

S  15  PPINT  "P-E  p-*"p:j;  "THETA. L-P  •  DEG-*"C 
SCO  W3*C 
S3 5  GOTO  'S15 
S 30  PPINT 

S  35  F* '  ■ A- 1  •  ■ A+l  •  •10.5 
340  G*FNH’  1  111  • 

345  h=U3-G+'?0 

350  M4»ATN‘  F  TAN', G  ■  ■ 

S' 5 5  U5* '  Uh  ' -  F-H 

sso  def  FNC-  X  '*'  ■  ■  A+i  ',-3  •♦xt3  •  ■  i  +  •  •  a- i  3  ■ 

SS5  1 0»FNC ' W 1  ■ 

370  .'  =  10 
■375  GO  SUB  3 SO 
3  30  0 1 =L 
335  GOTO  S05 

s so  »  »>  <  i-::».  ■  <:-t- Ft2>  >to.5 

S'? 5  L=ATN'  (.  >.-  F-ATN'  K  F  > 

SOO  RETURN 

SO 5  U*=-3-F  •  •  C OS •.  W4  •  t  *,  <  t.  3^A»-1  •  •  3-'  A-  1  •  •  -  -  3  IN'  W4  •  rO.  5  • 

SIO  GOTO  1030 

S15  IF  DT="YES"  THEN  SSO 

S30  FF'INT  "NOZZLE  IS  NOT  IDEAL  CONICAL" 

S3 5  DISP  "ENTER  NOZZLE  LIP  MACH  #='! 

SSO  INPUT  B 
S  35  U l =6 

S40  FPINT  "NOZZLE  LIF'  MACH  NUMBER* 'E' 

S45  DISP  "ENTEP  EXIT  VALUE  OF  U1-*'' : 

S50  INPUT  CO 
S55  U3 *C  0 

SSO  PPINT  EXIT  VALUE  OF  Ul*  'FROM  M. 0. C .  ■= "C 0 
SS5  DISP  "ENTER  RADIUS  OF  IMPINGEMENT.  WALL*"? 

S70  INPUT  0 

S75  PRINT  "RADIUS  RATIO.  IMPINGEMENT  WALL*"® 

S7S  Ct 101*0 

see  goto  S38 

935  GOTO  1030 

990  PPINT  "IDEAL  CONICAL  SOURCE  FLOW  NOZZLE" 

995  e*Wl*FN0<R8*RS*2  l  *C0$(W3)  >  > 

1000  PRINT  "NOZZLE  EXIT  MACH  NUMBER*''Wl 
1005  C0*W2*3*SIN<W3)*SQR<FMC<W1 >>/<Wt*WJ) 

1010  DISP  "ENTER  RADIUS  OF  IMPINGEMENT.  WALL*" . 

1015  INPUT  Q 

1020  PPINT  "RADIUS  RATIO.  IMPINGEMENT  WALL*'P 

1022  CC I01*Q 

1025  GOTO  830 

1030  W7»W2*H6 

1035  CC  I  >A 

1045  CC  2  3*W1 

1055  CC  3  )*C 


ji 


✓ 


£ 


c~ 

« 


r* 
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J 


£ 


« 
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1066  DISP  'ENTER  RE-DELTA1*""  5 

1067  INPUT  C£3J 

1068  PPINT  "PE-DELTfll“"Ct  5  ] 

1069  DISP  "ENTER  DELTA1- P*=" « 

1070  INPUT  D 

1071  PPINT  "DELTA1/R** "D 

1072  DISP  "ENTER  8.L.  EXP.  N="! 

1073  INPUT  N 

1074  CE15]*N 

1073  DISP  "ENTER  FIRST  ESTIMATE  JO*"! 

1083  INPUT  JO 

1093  PRINT  "FIRST  GUESS.  JO-FACTOR*" JO 

1103  cmi*  jo 

1110  GOTO  1165 
1115  I'EF  FNO<::> 

1  ISO  M 1  *  1 . 4 
1125  N2*FN J1 Ml ' 

1130  IF  ABS(:;-M£'-  0.0001  THEN  1155 
1135  M4* ■  A*  1  ••Ml  • M 1 -M 1 ♦ • A- l  ■  *2 ,•  -  1  Ml 
1140  M3*M1*»  X-M2  1  •  X.-M4  • 

1145  M1*M3 
1150  GOTO  1125 
1155  P ETUPN  M3 

1160  DEF  FNJ*  ::>*>,  v2  ■  A+ 1  ■  1  +  ■  A-  1  :  2  ■  •  T  •  •  A+ 1  >  •  2*.  A-l 

1165  CC 12  ]=P3 

1167  PRINT  "  B .  L .  E::p.  N*"N 

1168  CC  6  ]*D  PS 
1170  LOAD  *5. 3 
1500  END 


I 

r.' 


1  \ 
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1  COM  CC  25  ] 

10  PPINT  JET  IMPINGEMENT  ON  CYLINDRICAL  HALL*  HHK ■  3  25  80 
15  DIM  AlC  3  ]•  ETC  3). CtC  3] 

20  DISP  DATA  ENTERED  6Y  DIGITIZER*""!" 

25  INPUT  El 

30  IF  61# " YES ""  THEN  290 

35  PPINT  "DIGITIZATION  OF  STRIP-CHART  DATA" 

40  PPINT 

43  PPINT  "DIGITIZER  COORDINATE  ORIENTATION" 

50  PPINT 
55  ENTER  ■  9. 

60  TO*ATN'  y  ::  > 

65  PPINT  "ANGLE  THETA* "TO- 180  PI "DEGREES" 

70  DIM  6C 100« 2  ]> AC  100. 2  ]< DC  100- 2  1 

75  DISP  "ENTEF  ORIGIN  COORDINATES  OF  GRAPH.  ;:S.YS»"! 

SO  INPUT  ::S> v? 

35  PPINT  "ORIGIN  OF  GRAPH  AT  LOCATION  ::*"::35  "Y*"Y3 
90  DISP  "ENTER  X9.  REF.  VALUE  X  AXIS*"* 

95  INPUT  :.9 

100  PRINT  "REF.  VALUE.  X  AXIS*"!:? 

103  WRITE  ■ 9. *  • 

110  ENTER  ■  9 •  ■*  ' X0 *  Y0 

115  X'9*X9-X8 

116  PRINT  X9 

120  DISP  "ENTER  Y9.  REF.  VALUE  V  AXIS*"! 

125  INPUT  Y9 
130  Y9*Y9-Y8 

135  PPINT  "REF.  VALUE  Yl,  Y  AXIS*"Y9 
140  WRITE  < 9 1 * ) 

145  ENTER  <  9 «  *  > X 1 »  Y 1 

150  DISP  "ENTER  »  OF  POINTS"! 

153  INPUT  N0 

160  PPINT  ""NUMBER  OF  POINTS*"N0 

163  FOR  1*1  TO  N0 

170  WRITE  <9.0 

175  ENTER  <9. *>DE I  *  1  ], D£I . 2 J 

180  WAIT  200 

183  NEXT  I 

190  FOR  1*1  TO  NO 

1*5  Pf  T .  t  i*nr  l .  i  i-rnc.Ti-unr  x .  ■>  tkto 


w 


200  BC  I.2  3=DC  I »  2  3*C0ST0-DC  I,  1  3*SINT0 

205  ME!  T  I 

210  FOP  1*1  TO  NO 

215  «c  1. 1 3*x9*et  I >  1 3-  iOR'  ::o*::o+'.o*vo 

220  PC  I.23*YS*BC  I  >  2  3  SORO!1  *X1  +  Y 1  *Y  1  *♦'■'3 
225  next  l 

2  SO  FOP  1  =  1  TO  HO 

235  PPINT  "  I  =  "  1 5  fit  I.  1  05  AC  1-21 
240  he::t  i 

245  DISP  SHftLL  CRTS  BE  STORED*"  > 

250  INPUT  ft* 

255  IF  ft*#" YES"  THEN  3S0 
2 BO  STORE  DftTft  #5. 4. ft 

2*5  PPINT  "DftTft  STORED  IN  FILE  4. ON  FLOPPY  DISC" 

2 TO  PRINT  "NOTE  DIM'  100-2  ■“ 

275  PPINT  "IF  POINTS  TO  BE  PLOTTED  CONT  345"! 

2S0  PRINT 
2S5  STOP 

2 SO  D I IP  POINTS  FROM  t  EYBOfiPD*" ! 

235  INPUT  C * 

•00  IF  c  *# " YES"  then  345 
305  DI3P  "ENTER  #  POINTS*" * 

310  INPUT  N 
315  FOP  1=1  TO  N 

320  DISP  "ENTER  "I"  •* Y-  "I"  ■  =  "! 

325  INPUT  ftC I < 1  ]>  ftC I  -  2  3 

330  ppint  •=  ■  i :  ••::=  pic  I <  i  is  "v="ft[  I •  2  J 

3  35  NE  T  I 
340  GOTO  245 

345  PRINT  "DftTft  LOftDED  FROM  Dili  FILE  #4" 

350  LORD  DftTft  #5- 4. ft 

355  DISP  "NUMBER  OF  POINTS  IN  FILE'="i 
3*0  INPUT  N 
3*5  GOTO  375 
370  N»NO 

375  FOP  1=1  TO  N 

3  SO  PPINT  "I  =  RC  1.1  35  ".  =  RC  I  o  2  3 

3S5  ne::T  I 

3 SO  PRINT  "SELECT  MftX  RND  MIN  VftLUES  SUCH  THRT  DIFFERENCES  RFE  FCUER 
395  PPINT  "OP  2  OP  5  TIMES  10TN" 

400  PRINT 

405  DISP  "ENTER  : :-M i n. : : - m ft : : . y-m i n. y-mr: : " : 

410  INPUT  Pi  <P2<01 <02 

415  SC  RLE  -0. 1 *P2  + 1 . 1*P1 •  1 . 1-P2-0. 1-F !• 1 . 1-01-0. 1 *02 • 1.1*02-0.1*01 

420  ::r::i l  oii'.PS-Pl-  io*Pi<P2 

425  FOR  : = P 1  TO  P3  STEP  'P2-P1'  10 

4  30  plot  -  0 « 1 

4  35  C PLOT  -4.-2 
440  LRBEL  ■ 4S5>X 

445  next 

450  Yft: :  I  '3  PI.  '02-01'  1 0.01  >  02 

455  FOP  .'=01  TO  02  STEP  >  02-01  >  10 

4*0  PLOT  P 1 • Y • 1 

4*5  CFLOT  -7. -0.3 

470  LRBEL  ' 485;Y 

475  NEXT  Y 

480  LRBEL  (*• 1.5. 1.7.0.8- 10> 

485  FOPMftT  F6.3 

490  DISP  "ENTER  #  POINT  PAIRS*"? 

495  INPUT  N 
500  FOR  1*1  TO  N 

505  PRINT  "I*"I» "X="ftt I» 1 35 "Y="ftC 1.2 3 
510  PLOT  ftC I. 1  1. At  1-21 
515  CPLOT  -0.3. -0.3 
520  LRBEL  <#)"0" 

525  IPLOT  0.0 

530  PEN 

535  NEXT  I 

540  STOP 

545  P*CC 103 

555  PPINT  'R,U/R.E*"R 


,v?'JV.w??y?>.9fi'!‘r>w.m?.9>:\,Ji.aP.m*7*r\wyr*rp'\w’ v-1* •y.n 


£r 

,*.V  . . 

,*V 
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5 <50  K-CC  1  3 
579  PRINT  "K*"K 

575  DEF  FMA<Y>*<<K+l>-'.k-l>*Y-*Y'>  ■  2-Y*V  • 

530  Y«CC16J 
530  m»Y 

595  PRINT  " M * ' 0 * " Y 
6 ft 0  R'?=FNP  ■  Y ,' 

605  B*CC 17  3 

615  PRINT  “DM*  DR*  REATTACHMENT="B 
61ft  A*l-FNAtY> 

625  PRINT  "Ul+*"-2*A*B 

6  s  0  L 

6'5  D I SR  "ENTER:  INTEGRATION  STEP  DELTA-7  R  •  u='  : 
640  INPUT  & 

►.45  FF I  NT  "INTEGRATION  STEP  DELTA  ::  R.U=D 

650  ::o*cus3 

660  plot  70*rs 

665  FOP  1=1  TO  to 

67ft  H2=m-'  C  •  2*>'  1  -FNA1  HI  '  '  >  '*D 

675  P2=FNP‘ M2  ’ 

6  a  ft  PRINT  "V  R,m«"!-R*D»"R  P0.M*"P2? “H*»  112 

625  ::=::0*I»P*D 

630  PLOT  ::*P2 

6S5  rn=M2 

700  ne::t  i 

705  REN 

7io  : : i =  c c 7 ] 

720  F  i=ar?3 

720  print  r.e.  '•  1 2C.  =  -::i R*E-INVI2C.  =  "::0 
7*5  PRINT  " P >  B  RO* N*  ”R  1 5  "R* F'EAk  PO*N*"P? 

740  DEF  FnRi  ::.'  =  ■  l-'  i -i  -k+1  >  *f*'f 

745  PLOT  0«P1 

750  plot  : :  1  *  F 1 

755  PLOT  70  •  F'S 

760  PEN 

765  END 


1  CON  CC253 

10  DIM  AC7, 753<BC2<  20  3<1  C2-753 

20  PRINT  "THIS  VERSION  WITH  PLOT*HH»  2  20  SO" 

20  PRINT  "LINEAR  MINING  PROFILE  IN  EXPANDED  I  N  PROFILE  S.L. 
40  PRINT  THIS  VERSION  OF  8.  23  80  CAN  HANDLE  THIN  ARPR.  B.L. 
50  DEF  FNZ  •:  F  >  =FTN  SOP <  1  -  ■  1  -C 1  T2*F  1 2 > *P  t  ■.  •  G - 1  •  G  •  • 

60  DEF  FNC  •  M  1  =SC!R 1  NT  2  v 2  • G- 1  ■  +MT 2  •  • 

70  DEF  FNM>.  F  '  *FTN*FNDF  -  '  1 -Ft2+C IT  2  ' 

SO  DEF  FND '.  F  '•  =  <  1  -  •  l  -C  1 T 2+FT2> *Pt «,  •  G- 1  •  G  •  •  1 0 . 5 
SO  DEF  FNA'  F.'*FtlWl-Ft2*Ctt2> 

IC'O  DEF  FNYCF  >*02*' -F.3+L0G<  C  1+C0  >  ‘1-C0*  •  '2-C5  •  • 

110  FLOAT  4 

120  G=CC  1  3 

125  N=CC 153 

120  M1*CC  2  3 

135  N2=Ct 14  3 

140  D=10 

145  Pl=CC 5  3 

150  I =CC  7  ]/Ct 6  3 

234  PRINT  "X/DELTA.  FINAL*" I 

240  K*12+2.76*M2 

250  20*YO*F8* J*1 

2 60  P9* 1 0t  <  4 . 754+0.  0304112+0 .  05 3»M2 T2  -> 

270  P*<<l  +  <G-1>*M1T2/2V<1  +  <G-1>+M2t2.  2>>T<G- • G- 1 ' > 

280  PRINT  "GAMMA*"GI "UPSTREAM  B.L.  PROFILE  EXP.*"N 
290  PRINT  "Ml* "Ml I "M2* “M2 5 "P2'Pl*"P 

300  PRINT  "REYNOLDS  # »  < 1 -A-DELTA1 >*"R1 5 "SIGMA  CCONST.  '»"k 
310  PRINT  "*  OF  INTEGRATION  STEPS*"DI" REYNOLDS  #*  TRANS. *"FS 
320  Cl»FNCMl 
330  C2-FNCM2 

340  J3»SQRCl-<l-C2t2>*EXR(L0G<(l+C2>  •  1-C2,>>  C2-2'-  C2 
350  PRINT  "SIMILARITY  VALUE  FOR  J0*"J3 
360  7i«P*Pt*et-ClT?,'tP'ri 


i 


J70  F«Y1»I1*J1*O*0 
3S0  XI *23*0 
J90  Z5*1E+10 
400  Fl*F 

410  IF  F1*F8  THEN  440 
420  h=d-  io 
4  JO  03*J 

440  fj=fi+i  •  ii 

445  G0JU8  450 

44c'  GO  SUB  940 

447  GOTO  410 

450  F9*‘  Fi+FJ  •  2 

430  IF  F9  1  THEM  10 JO 

470  h=FMBF? 

4  J0  IJ=I 1+FNHF9-* F2-F1  • 

4*00  Y2*Y  1 +P  t  ■  - 1  G.**M*C  t*FMZF?4*  F2-F1  • 

500  YJ*Y2 

510  Dl»ll*Pt' -1  '-M-C  1  * -1+LuG*  ■  1+FMBFl  •  •  1-FMDF1  •  ■  •Z-FtU-Fl  •  • 

5  JO  I'J=I  J»P  T*  -  1  '  *M**C  1  -  •  —  1  +LOG*  <  1  +FHDFJ  >  '  1-FMHFJ  ■  •  ■ 2-FHDFJ • 

5 JO  IF  t>l  =  0  THEM  550 

540  J7=J-Ff<-l  >+Cl+FMIiFWI  B1 
550  I=2»P T  —  1  •*C1*FNDF2*M  B2 
5*0  JJsJl+FMftF?**  FJ-Fl  • 

570  IF  D 1  =0  THEM  530 

530  Z3*J1*Z7 

5 JO  IF  B1=0  THEM  330 

300  OOs'jOF* ■  1  - ■  1  -  >.  FMIiF  1  •  1 2  ■  *£: '.F'JS  1 

310  i 0=00  C2 

* JO  ->2=  J1+FNHF3*"  F2-F  1 

*J0  0 ?*FMBF  1 

340  F 4 = Q 0  03 

350  Z9*Dt+'  -F4+L0G1  1  1+00  *  *  l-00>  ’  '  J-OJ  •  * 

330  1 0=  JJ- J 

370  IF  B1#0  THEM  330 

330  Z3=0 

330  CO® 3 OF'1  1-*  l-'  FHDFJ  •  TJ  '+E::FIO  > 

700  JO=C0  C J 

710  IF  YO* 1  THEM  730 

7 JO  FFIMT 

7  JO  FFIMT  "  PH  I  - 1  i.  I )  * "  F 1  ?  "  PH  I  - 1  ’  I  ♦  1  * =  "  F  J 

740  FFIMT  "PHI-J'  I  "FMEF  1 ,  C  J  5  "PH  I  -  J*  I  +  1  '  =  "FMDFJ  CJ 

750  FFIMT  " V2  V1«”Y2* "BELTp.  BELTfil  =  "BJ5  D1 

730  FFIMT  "JO*’JO?"YL  r'ELTfll  =  "YJ-BJ* 

770  FFIMT  "  JO-FHTIO*'' J0/J3 

771  FFIMT  1 1 ;  us  j i s  jj;fi;fj;::i;::j 

7.30  C  l  *FNDF2  . 

730  IF  J0=0  THEM  1450 
300  J3=C  J-'  1- J0T2+CJT2  I'J 
310  FS*C0- C3 
3 JO  Z2*'  Z5+Z6  >  ''2 

3  JO  Z4=IiJ-'  -FJ+LOG'!  f  1  +C0>  *  •  1-C0>  >/<2*C3>  > 

340  T  l=3C'FP I  +  •: 1  * C2T2- 1  >  * >■  Z4-Z'? >  •' <.  D2-D 1 ) 

350  E0*T1*< l-C012>*C2/<fl*(l-C2*C2>> 

851  IF  Z0=1  THEN  370 

830  FFIMT  "SHEAR  STRESS  FUMCTION*"Tl 5 “EDDY-VISCOSITY  FUNCTION* "E 

370  X2*X1+<Z4-Z3>*Z1'Z2 

880  GOSUB  1790 

890  IF  Y0*1  THEN  910 

900  PRINT  "X'DELTflt •  LPM.*"X2 

910  R2*<C2»Rl*X2*P*ai-Clt2>''<l-C2t2>>t2>/Cl 

920  PRINT  "REYNOLDS  *  R/-2.X  =  "R2 

930  IF  R2>R9  THEN  1330 

935  RETURN 

940  23*24 

950  25*26 

960  XI *X2 

970  11*12 

980  J1 *J2 

990  Yi*Y2 

1000  F1-F2 

1010  RETURN 

1020  Y3*Y2 


105 


1050  PRINT 
1040  X1*X2 

1O50  PRINT  "EDGE  OF  SHEAR  LAVER  REACHED" 

1052  G*J 
1060  PRINT 

1070  SO*-1*LOG<  <  1*C2>'  <  1-C2>  >  1  2*C2.‘ 

1O30  Si*H*Cl*I2-  P-V2*C2T2  '  l-CiTi:> 

1 050  Tl*Cl+N*J2'P-Y2+C2.'U-C2t'2> 

1100  S2»Y  3»C2 1 2/ «.  I  -C2t 2  ) 

1110  D2*‘Sl+S2>' SO 

1120  PRINT  "V2/DELTA1*“Y35 "DELTA  DELTA 1*"D2 
1130  T2*Y3*C2-  ■.  1-C2T2.' 

1140  L=2*C2*1  Tl+T2>-'D2 

1150  JO-SOPCI--  t-C2T2,-*E::F'L  ■  C2 

1160  PRINT  "JO«"J0 

1170  PRINT  " JO-RATIO*” JO  J3 

1130  IF  P2.R9  THEN  1630 

1130  26=C2*< t- JO* J0+C2+C2  •  D2 

1200  F2*C0  C2 

1210  22*' 25*26  •- 2 

1220  24=D2*( -F2+L0G'  • 1 +C01  • 1 -CO  1  •  '  2-C 2  •  • 

12  20  X2=::i  +  (  24-23  '*21  22 

1240  R2®1  C2*P1  *X2*F'  +  <  <  1  -Cl  t2>  '  •  1-C2T2  >  ’  1 2  •  Cl 
1250  PRINT  DELTA  I  •  LAM.  =  "X2S  " REYNOLDS#  P-2..  :  =  "P2 
1260  PRINT 
1270  :.i»::2 
1230  Y  2  =  Y  2  +  2 

12'?0  IF  P2  2E  +  06  THEN  1320 
1  2  0  0  G  0  3  U  E'  1 7  3  0 
1310  GOTO  1100 
1220  PRINT 

1320  IF  YO*0  THEN  1400 

1340  F3=F2 

1250  D*  1 0*D 

1260  06 *J 

1 270  Y0*0 

1330  20*1 

1330  GOTO  410 

1400  PRINT  'TRANSITION  HAS  OC CURED" 

1410  20*0*0 

1420  os*:  :2 

14 20  07* J 

1440  xi*x2 

1450  GOSUB  1550 

1460  PRINT  DELTA  1  <  TUPB.*"X2 

1470  Tl  =  :OPPI*<  1  C2T2-1  >*'  24-23  •  D2-D1  > 

1430  E0 = T 1  * <  1  rC  0 1 2 ,•  »C 2  ' C  A* C  1 -C2+C2  •  • 

1430  PRINT  "SHEAR  STRESS  FIJNCTION*"TU  "EDDY-VISCOSITY  Fi.INCTJi 

1431  PRINT  Ill  I2S  Jl?  J2?F1!F25X15!<2 
1500  GOSUB  17?0 

1505  IF  0*0  THEN  1550 
1510  GOSUB  940 
1512  GOTO  410 
1520  GOSUB  1550 
1530  GOSUB  1790 
1540  GOTO  1460 

1550  F3*C0-'C3 

1551  IF  0*1  THEN  1560 

1552  Z3*U4 

1553  GOTO  1570 
1560  23*24 
1570  24*FNYF3 

1580  >!2*X1*<24-23)*K*D2/CFNYF3*SGRPI) 

1581  PRINT  X1JX2.F1IF2 

1582  PRINT  "23*"Z3I  “24*"24! "F3*"F3! "CO*"C0J "C3='T3 

1590  RETURN 

1600  S2*Y3#C2t2/'<  1-C2T2) 

1610  D2-(S1+S2>/S0 

1620  PRINT  "Y2/DELTA1*"Y3? "DELTA/DELTAt*"D2 
1630  T2-Y3*C2'<l-C2t2> 

1640  L*2*C2*<T1*T2)/D2 
1 650  J0«SQR <1  - <  1  -C2t2 )  *EXPL  >  -'C2 


106 


1660  PRINT  •,J0«"J0 

1670  PPJNT  -  JO-RATIO*'^©  J3 

1600  F 3*CO'C2 

1600  23=24 

irOO  24=FNVF3 

ir  10  : :2=x  1  +  < 24-2 3 .■  +i  ,  -D2  < fmvf 3+sopp i  > 
in 3  K  1.  Jl*:'2 

1713  1C2«  J3*J0 

1720  F'2*'  C2+P1-X2-P*1  •  1-C1T2'  "■  1-C212.1  >12'  Cl 

1730  PPIMT  ''X' DELTA  l  >  TUPE.  =  "X25  “REYNOLDS#  P  2-X='P 

1735  IF  ::2M  them  2065 

1740  GOSUB  1700 

1750  V 3='i' 3+0. 5 

1760  xi*::2 

1770  IF  P2  2E+07  THEM  1020 

1750  GOTO  1600 

170O  K 1. J3=AC l.JJ*X2 

1300  AC2.JI-D2 

1310  Iti-  Jl-AC  3*  J  3=  JO 

1320  fit  4. J  1*1 *X2  '• E0  +  D2 > 

1330  AC 5* J 3*24 

1331  AC  6*  J  ]*I2 
13  32  AC  ?«  J  3=  J2 
1340  FEINT 

1350  PPIMT  " J*" J 

1352  IF  :2  I  THEN  2710 

1360  PPIMT 

1330  00= J 

1S0O  J-.J+1 

1000  F’ETUFM 

1010  GOTO  1600 

1020  : TOP 

1030  P=- 1  NT 

1040  FORMAT  3X.  "  J"  >  3X.  ''X  DELTA.  1  2::- "DELTA  delta.  1 

1050  WRITE  ‘ 15.1040' 

1060  FOP  J-l  TO  00 

1061  IF  J  =  07  THEM  1070 

1062  AC  4»  J  1*0 

1070  UFITE  '.15.  1000  'J.  AC  1  •  J  I .  AC  2  •  J  ] « AC  '3  •  J  I  *  At  4  <  J  1 
1050  XT  J 

1000  FORMAT  1F4.0.4E12.2 

1  005  FPInt  '  NOTE  SCALE  LGT .61* LGT5000* 2>  2" 

2000  GOTO  2210 

2010  FOP  J=07  TO  00 

2020  plot  lgt-  ac  i«  ji.-.ac  3.  jj 

2030  NEXT  J 

20  35  LABEL  "  JO" 

2040  FEN 
2050  5 TOP 

2060  FOP  J=07  TO  00 

2070  PLOT  LGT • AC  1  *  J  3  >  *  AC  4  * J  I  1 0000 

2030  NEXT  J 

2085  LABEL  •.*>  "PE-K.  T/10E4" 

2000  PEN 
2100  'STOP 

2110  FOR  J*Q7  TO  Q0 

2128  PLOT  LGT<AC1»J3)«AC5.J]  20 

2130  NEXT  J 

2135  LABEL  <*>'Y2'20" 

2140  PEN 
2150  STOP 

2160  FOR  J«Q7  TO  <39 

2170  PLOT  LGT<AC1» J  1>»  At  2*  J  3'20 

2180  NEXT  J 

2185  LABEL  <*> "DELTA/20" 

2100  PEN 
2200  GOTO  2670 

2210  DISP  "ENTER  LOWER  LEFT*"* 

2220  STOP 

2230  DISP  "ENTER  UPPER  RIGHT-"! 

2240  STOP 

27*0  ®rftl  F  I  r,TA.  01.I  GT5A0A. -«.?•> 


r» .n 
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2260  :-:«XIS  0-LGTIO.LGT0. 1-LGT1O00 

2270  '.'AXIS  LGT0. 1.0. 1.0. 1.3 

2230  A*0. 1 

22*0  PLOT  LGTA.0 

2300  PLOT  LGTA. 1.5 

2310  FEN 

2320  0*10*0 

2330  IF  A  10CO  THEN  2250 
2240  GOTO  22*0 

2  250  FOR-  E *0  TO  1.5  STEP  0.5 

22-50  PLOT  LCT0.1-B 

2270  PLOT  LGT 10O0. 8 

2330  FEN 

23*0  NEXT  6 

240O  M=0. 5 

2410  PLOT  LGThi.0 

2420  plot  lgta-i.5 
2420  FEN 
2440  A* 1 0-A 

2450  IF  A  10OO  THEN  2470 
2430  GOTO  2410 
2470  : TOP 
2430  *  T ANDAP  D 

24*0  fop  x»-i  to  2 
2500  : := inT: : 

2510  plot  -o.  1 . 1 

2520  CPLOT  -2  <  0 

2520  LABEL  •  *  ■ "  10 1  “  InT:: 

2540  next  :: 

2550  FOP  1*0  TO  1.5  3TEP  0.5 
2550  PLOT  -1.5.Y.1 
2570  CPLOT  -0.2 
2530  LABEL  ' 

25*0  NEXT  y 
2500  FORMAT  F3.2 
2510  PLOT  -0.5- 1. 5* 1 

2520  LABEL  '-'"FREE  3H6AP  LAYER  AFTER  Ex  L.  EXFAHi I 3n“ 
2520  LABEL  <>•  “BOUNDARY  LAYER.  £XP.“N'.  RE. DELTA*  1  *  F'l 

2540  LABEL  •'* '"APPROACH  MACH  ••“Mr  JET  MACH  *»"M2 

2541  LABEL  i*  ••,GAMMA«,  G"?IGt'1A*  "1 

2542  PLOT  3.05-0'  1 

2544  LABEL  •*•''::  DELTA-  1' 

2550  GOTO  2010 
2550  END 

2570  PLOT  LGT' AC  1-07  ]  -  -0 
2530  PLOT  LGT' AC  1 • 07  3 ' » 1 . 15 
2530  PEN 

27C0  LABEL  <*•  1.5-  1.7*  PI  2. 3  10,'  “TRANS I  T  ION” 

2701  PLOT  LGT ( I > . O 

2702  plot  lgtci>.0.35 

2703  FEN 

2704  LABEL  C*  •“RECOMPPESSIOH  START" 

2705  GOTO  2020 

2710  IF  0*0  THEN  2820 

2715  IF  F 1  >0.3  THEN  2355 

2720  W1 «X1 

2730  U2-X2 

2740  V2«F2 

2750  Vl-Fl 

2760  Ul-Il 

2770  U2-Y1 

2780  U3-J1 

2785  U4»A[ 5»  J*1 3 

2790  W0-AC 1. J3-I 

2795  U5-C3 

2000  H»0“1 

2810  Bt  2< 1 3*At 1 » J-l 3-1 
2815  Bt 1.1 3«R1 

2820  BC  1-HM  3«V2-W0*<BC  l.H3-V2)/<BC2.H3-W0) 

2822  PRINT  “H«“MIBt 1  *  H+ 1 31 BC l >  H  3! BC  2- H  3 
2830  Fl-Vl 
2840  ll-Ul 


7 1  -U2 
J1-U3 
>u»ui 
C3-U3 

F2-BC l.H+1  ] 

COSUB  450 

IF  ABSF  Bt  2 >  H  I  1  IE-0  3  THEN  2S5U 
H*H»  1 

bc2»hi«x|-i 

p'pINT~"S0LUT10N.  X2="Bt  2*H>T*  "  JO-“  JO-* "  JO-RAT 10-"  JO 
IF  ABS'CUU-JO -  J3  •  0.01  THEM  29ftU 
CC  1 1  3= JO  J3 
LOAD  »3*3 

fpjilr  •'THIS  IS  THE  FI  URL  RESULT  FOP  THE  SUB  I1  EM  ENLARGEMENT  F  f  .  EtEM 

'jO-KS*  J-1  ]*‘H:  2.  J3-KC2*  J-l  3  I-KC  I.  .1-1  3  •  •  1 1  1  •  J  3-fCl  •  1  j  • 

i  PRINT  "JO*"JO“AT  :xBELTAl-"I 

;  cc  n  i= jo  J3 

:  GOTO  2S52 

r  pi SR  "IF  plotting  desired,  enter  0“« 

I  INPUT  U? 

\  IP  U?#0  THEN  SOSO 
!  GOTO  IS  SO 

)  lord  #3* I 

>  END 


COM  CC  25  3 

PRINT  "PLUME  IMPACT  ON  CYLINDRICAL  MALL*  2EP0  6L0M-6Y.  THEORY  AND  E 
3  PRINT  "PLUME  IMPINGEMENT  ON  R  CYLINDRICAL  MALL*  HHt  . i  21  7S" 

J  PRINT  "ITERATIVE  SOLUTION  FOR  RINDING  ZERO  BLOU-BY  CASE*  HHt  S  5  “V 
3  PRINT  "CALCULATION  OF  APPROXIMATE  FLUME  SHAPE" 

5  PRINT  "FOR  IDEAL  OP  NON-IDEAL  CONICAL  C-B  NOZZLES" 

3  DIM  DPI 3  3*113*33 
5  DEG 

3  DEF  FNA>  ',;  ■  aRTNO.  SOR'  1  1 E-S  ? .•  • 

5  DEF  FN6  <.  X  >  «ATH<  SOR  <.  1  >X  >  C+1E-S?  • '♦2-ATNlE+S?-' X  0* 

3  DEF  FNP<M  >*<• 2*< A-t >-M-*M>  ■ A+I1  ■  (••  ■ A*1  ■  <***A-l  •  •  •  M 
5  Z=I 
3  "=l 

3  PRINT  "NOZZLE  FLOW  SPECIFIED" 

3  A-CI I  3 

3  PRINT  "GAMMA- "A 
5  M 1 -B-CC 2  3 
3  C-CC S3 
5  CO-Ct 4  3 
3  F3-CC 12  3 
3  M3-C 
3  GOTO  333 
3  PRINT 

3  Pa. <A-l>  ' A+ 1 >  > TO. 3 
3  G=FNA< I /U 1 > 

5  H«U3-G*90 
3  W4«ATN<F/TAN<G>> 

5  U5»<W4VF-H 

3  DEF  FNC(X)-<  <  <  A+ 1 2>*Xt23/,(  1  ♦(  <  A-l  >*.X»X)/'23 
3  I0«FNC<WI> 

3  K-10 
5  GOSUB  890 
3  01  «L 
3  GOTO  903 

)  K-<<1-X>/(X-1/FT2))T0.3 
3  L-ATN(K)/F-ATN<K/F> 

3  RETURN 

5  W6a-2*F,<C0S<W4>r<<<3*A)-I V<2#<A-l>>>*SIN<W4>t0.5> 

3  GOTO  1030 
3  W1«B 

3  PRINT  "NOZZLE  HP  MACH  NUMBER* "B 
3  CO-Ct 4 3 
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955  W2*C0 
970  C!*CC  101 

975  PRINT  "RADIUS  RATIO.  IMPINGEMENT  UALL*"Q 
988  GOTO  830 
935  GOTO  1030 

995  B*W1=FN0(R8*R8*2''<  1+C0S<W3>)> 

1OO0  PRINT  "NOZZLE  EXIT  MACH  NUM6ER="H1 
10O5  C0*W2*3+SIN<M3>*S0R<FNC'm  >.'•■  (W1*W1  • 

1O10  D I SP  'ENTER  RADIUS  OF  IMPINGEMENT.  WALL="S 
1015  0=CC 103 

1030  PRINT  "RADIUS  RATIO*  IMPINGEMENT  WALL*"0 
1035  GOTO  330 
1O30  U7*W2*W6 
1O40  D=  1 . 3-W 1 

1045  PRINT  "FIRST  TRIAL  JET  MACH  NUMBER- 

1055  J0*CCm 

1OS0  PRINT  " JO-FACTOR® " JO 

1035  IC1.1]=D 

1O70  PRINT  "JET  SURFACE  MACH  NUMBER* “D 
1075  J=FNC  (  D.* 

1030  ::=  j 
1035  GO SUB  390 
1O90  03=L 

1095  E2*U4+F+<02-01-FNAa  D.+FNA'  l  6)  • 

110O  Dl  =  ’ E2-W4)  '10 
1105  B0=U5 
1110  El =H4 
1115  C 1 *W7 
1130  GOSUB  1130 
1135  GOTO  1190 
1130  SI =33=33=34=0 
1135  FOR  N=1  TO  10 
1140  E3*El+Dl*(N-l/2> 

1145  S1=C0S<E3 'FUSIN' E3' t  ■ -0.5 '-COS' ES  't  '  -l  ^  2-F T2  ■  '-D1  + 
1150  S3=SIN<E3-'F'.*SIN'  E3 'T'.:-0.5.'-C0S' E3:.'t' -1  »2*Ft2'  '*D1  + 
1 1 55  S  3  *  3 1 N '  E  3  /  F  >  *  S I N  ■  E  3  :•  t '  - 1 . 5  '  * C  0 S  *  E  3  >  T  (  ...  A  -  3  •  • 

1130  S4=C0S'.E3-  F)+$IN<E3>  tv-1.  5>+C0S'.E3  'Tv  (A-3)  •  3*'  9-1  .  . 
1135  NEXT  N 

1170  U O * — ■  •  C0S<E2>t<  < 3+A-l  •.  <  2*' A-l  >  >  >  >»SIN<E2)t0.5  •  •  3*F 
1175  U1*U0+'  •  S 1 - F 4 S 3 > *  C  0  S  <  B  O  >  + ' .  S  2 + F  ♦  S  4 > -SIN' B0'+C 1-130  PI 
1130  PRINT  " IJ 1  * 1  F > * " U 1 
1135  RETURN 
1190  G2*FNAC  l  D  ;■ 

1195  T3=C +02-01 

1300  R0=-<  •  U 1  -  SO R <  J > + S I N •  T 2 )  +S I tH G3 >12)  '.SIN'  3+G3  •  •  • 

1305  Pl*-1  R8 
1310  PRINT 

1315  PRINT  "SOLUTION- 

1330  PRINT  "INITIAL  SLOPE  OF  PLUME.  THETA-F="T3 
1335  PRINT  "'INITIAL)  RADIUS  OF  CURVATURE.  PLUME*  PI 

1330  IF  R 1  * ( 1 -COS ( T2 ) )  + 1 >0  THEN  1235 

1331  D=D+0.2 

1232  PRINT  "NO  PLUME  IMPINGEMENT" 

1333  GOTO  1035 
1235  PRINT 

1240  PRINT  "PLUME  SHAPE  APPROXIMATED  BV  CIRCULAR  APC " 

1245  WRITE  (15.1250) 

1250  FORMAT  5X> "X(-)".3X» "R<-)“»5X. "THETR(DEG)  " 

1255  FOR  $*8  TO  1.2  STEP  0.1 
1260  T*FNA(SIN(T2)-(S'Rl>) 

1265  S*R1#(SIN(T2)-SIN(T) ) 

1270  R* 1 +R l * < COS ( T ) -COS < T2 ) ) 

1275  WRITE  (15. 1280>S.R.T 

1280  FORMAT  3F10.3 

1285  NEXT  S 

1290  PRINT 

1295  PRINT 

1300  G*FNA(  l/'D) 

1305  C*FN8(C0S(T2)+(Q-1)/R1) 

1310  B*D 

1315  E1*ATN(F/TAN(G)) 


pH 

110 

2 

J 

.  ,\ 

b‘,% 

»  . 

,  \ 

.‘In 

O  •. 

1325 

GOSUB  090 

1330 

02-L-C 

It  -  WBSSSm 

1335 

X4»FNC<B> 

_  _ 

1340 

X«X4 

j 

v> 

134? 

COSUB  090 

>  .v 

1350 

Ll-L 

** 

*.'> 

•  « 

1355 

::*X2*3.5 

•  «* 

L"  V 

1  3(50 

GOSUB  390 

A 

1  3(55 

L2-L 

> 

1370 

: : 3*’  02-Ll  •*•  ::2-::4  •  ■  L2-Li.'*::4 

•  . 

i  ri^a 

1375 

*  ‘  *  x  3 

1  .zM 

1380 

GOSUB  890 

s 

1335 

L3*L 

F^f 

S-; 

1  3?0 

IF  0B8CX3-X2  •  0.001  THEM  1420 

«* 

• 

,V 

1  335 

::4*::2 

1 400 

::2*::3 

1405 

L1-L2 

1410 

L2-L  3 

| 

1415 

GOTO  1370 

1420 

FF  IMT 

1425 

D«S0P*  2*X3/‘  0+1  -■  0-1  •*::3  >> 

;V 

14  30 

IF  :  o  THEM  1585 

■\ 

14  35 

MS *3 OF'  FMC  <  Ii  ,*  } 

1437 

CC 18 3-MS 

1440 

FFIMT  "M.  0TT.«"B*  "M*.  0TT . * "  $0R  <  FMC  ■  D  > 

,■/ 

1445 

C0»S0P..FMCvB>  •■•SIM.  G  2-I30S1  G-40  F1-8IM.  G '-SIM-C  <  • 

ul. 

1450 

C2=-2*F  -  C0S<EI  :•»  <...■  3*0  -  l  v  •<  2*' 0-1  ,.  >,+SIN(£l  •  r0.5  • 

• 

1455 

C 1-C0-C2 

& 

14(50 

BO-OTNCF  T0N(FM0'  1.  B  •  >  .•  F-90+FM0'.  1  •  B,' 

.,4 

( 

14.55 

E2*F*. 90-FMOc 1  D%+B0> 

-.**• 

1470 

t'l»- E2-E1  10 

•  " 

1475 

GOSUB  1130 

.  •*. 

■ 

14:30 

08«M8*U1  ■ 2*C0S> FMO1  1  D  •  ■  f  2 ■ 

kv. 

1405 

FPIMT 

‘A 

1  490 

PPIMT  "I  MV  ISC  ID  FEOTTOCHMEMT.  UOLL  PRESSURES" 

1 

1495 

P3».  1-. 0-t ‘*M8*MS  • 0* 1 ■  •  t • 0  '0-l  *> 

2  500 

PRINT  “  DM  -  -  DP  * "  OS  f "  F'EOK  PPESSUPE  POTIO-'FS 

fc*i( 

l.y 

1504 

CC  20  3-P8 

PV  ^ 

1505 

04-F 1 *( SIM'  T2 '-SIM'  C • • 

F\J 

150(5 

CC 1 7  3-Q3 

tr%* 

1510 

FFIMT  "L0C0TI0M  OF  PPESSUPE  FE01 •  ::.P  P-M02rLE="04 

•  «* 

K; 

1512 

1515 

CC 183-04 

FPIMT  "I MV I SC  ID  IMPINGEMENT  hNCLE«"C 

# 

H 

1520 

F'?«-08*P8*2*0*M8  •  C>»  -  -  0+1  ■-•  0-1  .*MS*M8  •  • 

1525 

FF  IMT  "INITIAL  SLOPE-  D •  P  PO )  D •  X-R - N :• « '' P9 

v 

1530 

C 3*0*0/. 2  -  ft- 1 > 40*0 > 

V 

1535 

Ct  -SC'P1  C  3  • 

ft 

1540 

E9«SC'P<Pt  1--  1  C  3-1  ■•  *•  -I+LOG'.  •  l*C6  •  •  t -C ■&  •  •  •  2-C6  •  •  ■ 

'  *- 

v 

& 

ft 

1545 

ES*<  •  T2-C -*PM30>*E?  '12*2.76*B> 

1550 

C3-FMO' SIMCC >+E8 '  SIM'C  V) 

1555 

15-50 

FPIMT 

Z*-I 

\> 

• 

15-55 

PPIMT 

1570 

FFIMT  "BOSE  PRESSURE.  VISCID  SOLUTION  CFOGE-  MODIFIED." 

*. 

1575 

C  «CS*0. 59 

1500 

GOTO  1315 

\ . 

. 

1585 

C4*D*D/'<2/<0-l  )*D*D> 

V‘- 

1590 

C5*SQR<  <C3-C4)/< 1-C4) )/C6 

s  \ 

1595 

C7*J0*$QR< l-< 1-C3>*EXP<L0G<< l*C6>'(l-C6>  >/C6-2>  >  'C6 

1600 

PRINT  "PHI »  D*“C5I "PHI i J*HC7 

•# 

1605 

D8*<T2-C8>*SORPI*C5/<<12*2.56*B>*180> 

A 

hr- 

1610 

C9*PHB< (C0S<C8)-D8)/< l*<T2-C8)*PI*E9/( 180*( 12*2. 56*B> > >  ■ 

••■;• 

1615 
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RETURN  M3 
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